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Abstract 

Let K C E, K' C E' be convex cones residing in finite-dimensional real vector spaces. An 
element y in the tensor product E <g> E' is K £g> if '-separable if it can be represented as finite sum 
y = x i ® where xi £ K and x\ G K' for all I. Let S(n), Tt(n), Q(n) be the spaces ofnxn 
real symmetric, complex hermitian and quaternionic hermitian matrices, respectively. Let further 
S+(n), H + (n), Q + (n) be the cones of positive semidefinite matrices in these spaces. If a matrix 
A G TL(mn) = TL{m) ® TL{n) is H+(m) (g> i/+(n)-separable, then it fulfills also the so-called PPT 
condition, i.e. it is positive semidefinite and has a positive semidefinite partial transpose. The same 
implication holds for matrices in the spaces <S(m) <8><S(n), TC(m) <S>S(n), and for m < 2 in the space 
Q{m) ®S(n). We provide a complete enumeration of all pairs (n,m) when the inverse implication 
is also true for each of the above spaces, i.e. the PPT condition is sufficient for separability. We 
also show that a matrix in Q(n) ® <S(2) is Q+(n) ® 5+(2)-separable if and only if it is positive 
semidefinite. 

Keywords: positive partial transpose, separability 
AMS Subject Classification: 15A48 

1 Introduction 

Let K 1 K' be regular convex cones (closed convex cones, containing no lines, with non-empty interior), 
residing in finite-dimensional real vector spaces E,E'. Then an element w € E ® E' of the tensor 
product space is called K ® K' -separable (or just separable, if it is clear which cones K, K' are meant), 
if it can be represented as a convex combination of product elements v <g> v' , where v £ K, v' £ K'. 
It is not hard to show that the set of separable elements is itself a regular convex cone. This cone is 
called the K (g> if'-separable cone. 

The notion of separability is intimately linked with the notion of positive maps [3J,[S]. Cones of 
positive maps appear frequently in applications [10], [6] and are dual to separable cones [5]. Separability 
itself plays an increasingly important role in quantum information theory [5]. 

A particularly important case in optimization and Mathematical Programming is when the cones 
K, K' are standard self-scaled cones such as Lorentz cones or cones of positive semidefinite (PSD) 
matrices [T]. Let S(n), rt(n), Q(n) be the spaces ofnxn real symmetric, complex hermitian and 
quaternionic hermitian matrices, respectively. Let further 5+(n), H + {n), Q + (n) be the cones of PSD 
matrices in these spaces. If the elements of a pair of matrix spaces commute, then the tensor product 
of these spaces can be represented by the Kronecker product space and is itself a subset of such a 
matrix space. If such a product matrix is separable with respect to two PSD matrix cones, then 
it is necessarily PSD itself. In the case of H + {m) (g> i? + (n)-separability, where n,m E N+, there 
exists another simple necessary condition for separability, the so-called PPT condition [S]. A matrix 
in Tt(mn) = H(m) (8 H(ri) fulfills the PPT condition if it is positive semidefinite and has a positive 
semidefinite partial transpose. 

In the spaces S(m) (8 S(n) and Ti(m) ® S(n) the PPT condition reduces just to positivity, i.e. 
inclusion in the cone S+ (mn) or H + (mn) , respectively. This is because the positivity property of real 
symmetric or complex hermitian matrices is preserved under transposition. However, the cone Q+(n) 
is invariant under transposition only for n < 2. Therefore for matrices in Q(m) ® S(n) the PPT 
condition is stronger than just positivity. Moreover, it follows that the PPT condition is necessary for 
Q+(m) (g> S+ (n)-separability only for m < 2, while positivity is necessary for arbitrary (n,m). 

The importance of the PPT condition and the positivity condition is based on the fact that these 
conditions are semidefinite representable (i.e. in the form of linear matrix inequalities) and hence easily 
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verifiable algorithmically, in contrast to separability. It is then important to know in which cases these 
conditions are actually equivalent to separability, rather than only necessary. In these cases one then 
obtains semidefinite descriptions of the corresponding separable cones. 

The theorem of Woronowicz-Peres states that in the case m = 2, n = 3 the PPT condition is 
not only necessary, but also sufficient for H + (m) ® H + (n)-separability [15], [8]. However, there exist 
matrices in 7i(2)<g>7i(4) which fulfill the PPT condition, but are not separable [TS]. Similarly, Terpstra 
[TU] has shown that positivity is sufficient for 5+(m) ® S'+(n)-separability for min(n, m) > 2, but not 
for n — m = 3. In [10] this was formulated in the equivalent form of sums of squares representability 
of biquadratic forms. 

One can then conclude that the positivity condition is equivalent to S+(m) <8> 5'+(n)-separability if 
and only if min(n, m) < 2 and that the PPT condition is equivalent to i?+(m) (g> -H+ (n)-separability if 
and only if min(n, ni) = 1 or m + n < 5. 

In this contribution we provide a similar classification for the spaces H(m)<g>S(n) and Q(m)®S(n). 
We show that positivity is equivalent to H + (m) (g) S+ (n)-separability if and only if m = 1 or n < 2 or 
m + n < 5, and that the PPT condition is equivalent to Q+ (2) <g> S+ (n)-separability if and only if n < 3. 
Further, we show that for m > 3 positivity is equivalent to Q+(m) ® 5+ (n)-separability if and only if 
ri < 2. In addition, we enumerate all pairs (n, m) for which the positivity property in Q(m) ® S(n) 
is preserved by the operation of matrix transposition, namely, the cases m = 1, n arbitrary, and the 
cases m = 2, n < 2. This involves mainly the following new and nontrivial results. 

First, we show that a matrix in Q(2) (8 5(3) is Q+(2) S'+(3)-separable if and only if it fulfills 
the PPT condition. Second, we provide an example of a matrix in H(2) ® 5(4) which fulfills the PPT 
condition, but is not H+{2) £g> S'+(4)-separable, thus sharpening the counterexample provided in [15] for 
the H+(2) (8> H+{A) case. Third, we show that if a matrix in Q(n) <g> S(2) is positive semidefinite, then 
it is Q+(n) C§> S'+(2)-separable. In addition, we provide examples of matrices in Q(3) and Q(2) <8> 5(3) 
which are PSD but whose transpose is not PSD. 

The remainder of the paper is structured as follows. In the next section we provide exact definitions 
of separability and of the PPT condition and consider some of their basic properties. In Section 3 we 
consider low-dimensional cases and relations between the cones we deal with. In the next two sections 
we prove the sufficiency of the PPT condition for separability in in the space Q(2) ® 5(3). In section 5 
we also provide an example of a PSD matrix in <2(2) <g) 5(3) whose transpose is not PSD. In Section 6 
we provide a counterexample against sufficiency of the PPT condition for H+(2) ® S+ (4)-separability. 
In Section 7 we prove the equivalence of positivity and separability in the space Q(n) ® 5(2). Finally 
we summarize our results in the last section. In the appendix we list facts about quaternions and 
quaternionic matrices which we use for the proof of the main results of the paper. There we provide 
also an example of a matrix in <5+(3) whose transpose is not PSD. 

2 Definitions and preliminaries 

In this section we introduce the cones we deal with and provide definitions of separability and the 
partial transpose. For basic information related to quaternions and quaternionic matrices we refer the 
reader to the appendix. Throughout the paper, i,j,k denote the imaginary units of the quaternions 
and the overbar 7 the complex or quaternion conjugate. 

Let further eo, . . . , e m -\ be the canonical basis vectors of R m . By ids denote the identity operator 
on the space E, by I n the n x n identity matrix, by nX m a zero matrix of size n x m, and by 0„ a 
zero matrix of size n x n. Let further diag(A, B) denote a block-diagonal matrix with blocks A and 
B. For a matrix A with real, complex or quaternionic entries, A* will denote the transpose, complex 
conjugate transpose or quaternionic conjugate transpose of A, respectively, and rkA the rank of A. 
Further we denote by GL n (R) the set of invertible matrices of size n x n with entries in the ring R. 

We now introduce several convex cones we deal with. 

Let L n be the n-dimensional standard Lorentz cone, or second order cone, 

i„ = |(x ,...,x„_i) T GR™|a;o> \[x\ + ■ ■ ■ + x 2 n _^ . 

Let 5(n) be the space of real symmetric n x n matrices and S+(n) the cone of positive semidefinite 
(PSD) matrices in 5(n); TL{n) the space of complex hermitian n x n matrices and H + (n) the cone 
of PSD matrices in TL(n); Q{n) the space of quaternionic hermitian n x n matrices and Q + (n) the 
cone of PSD matrices in Q(n); Q k (n) the space of quaternionic hermitian n x n matrices with zero 
fc-component and Q\(n) the intersection of the PSD cone Q+(n) with the space Q k (n). All these cones 
are regular, i.e. closed, containing no line, and with nonempty interior. 
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Let E, E' be real vector spaces and A C E, A' C E' regular convex cones in these spaces. 

Definition 2.1. [5] An element of the tensor product E ® E' is called A £g> K' -separable if it can be 
written as a sum J2i=iPi x i ® where A € N and p/ > 0, S A, xj S A' for alH = 1, . . . , JV. 
The set of A ® A'-separable elements forms a regular convex cone in E ® E' . 

We are interested in the case when the spaces £7, £7' are spaces of hermitian matrices, and the cones 
A, A' are the corresponding cones of PSD matrices in these spaces. If the elements in the factor spaces 
commute, then we can represent tensor products of matrices by Kronecker products, which will again 
be hermitian matrices. If the elements do not commute, then the Kronecker products will in general 
not be hermitian. Since we are interested in describing the separable cone by a matrix inequality, we 
do not consider this latter case in this contribution. 

In particular, for m,n € N+ we consider the space S(m) <g> S(n) as a subspace of S(mn), the spaces 
H(m) (g) S(n) and H(m) <g> H(n) as subspaces of H(mn) and the space Q(m) <g> S(n) as subspace of 
Q{mn). As can easily be seen, exchanging the factors in the tensor product is equivalent to applying a 
certain permutation of rows and columns to the corresponding Kronecker product. Hence exchanging 
the factors in the tensor product leads to a canonically isomorphic space, and we can restrict our 
consideration to the cases listed above. 

Definition 2.2. An automorphism A of E is called an automorphism of A if A[K] = A. 

The automorphisms of a cone A form a group, which will be called Aut(A). For regular convex 
cones A, A' let g £ Aut(A),<?' € Aut(K') be elements of their automorphism groups. Since g,g' are 
linear automorphisms of the underlying spaces E, E', we can consider their tensor product g®g', which 
will be a linear automorphism of the space E (g> E' . The following assertion is trivial but nevertheless 
very useful. 

Proposition 2.3. Let Qk®k' be the set of maps {g ® g' : E ® E' — > E ® E' \ g £ Aut(K),g' £ 
Aut(A')}. Then Qk®k' is a group. It is canonically isomorphic to the factor group [Aut(A) x 
Aut(A')]/{(a idy , a -1 idy) | a £ The cone of K K' -separable elements is invariant with respect 

to the action of Qk®K' , and Qk®K' is a subgroup of the automorphism group of this cone. □ 

Definition 2.4. A face F of a convex cone is a subset of A with the following property. If x, y £ A 
and £ F, then x,y £ F. For x £ K, the face of a; in A is the minimal face of K containing x. 

It is not hard to see that a face F is the face of a point x if and only if x is contained in the relative 
interior of F. If x £ K and y £ A', then the face of (x, y) in K x K' C E x E' is F x x F y , where F x is 
the face of a; in AT and F y is the face of y in if'. 

Definition 2.5. An extreme ray of a regular convex cone A is a 1-dimensional face. Any non-zero 
point on an extreme ray is a generator of that extreme ray. 
A convex cone is the convex hull of its extreme rays. 

Lemma 2.6. Let K C K. be a regular convex cone and L C M. N a linear subspace of dimension n. 
Let K' = K n L and let x be the generator of an extreme ray of A'. Then the face of x in K has at 
most dimension N — n + 1. 

Proof. Let d be the dimension of the face A of £ in A. Then the intersection F n L has dimension 
6! > d + n — N and x is in the relative interior of F n L C A'. On the other hand, d' = 1, since x 
generates an extreme ray of A'. It follows that d < N — n + 1. □ 

Any face of S+(n), H + (n), or Q+(n) is isomorphic to the cone S+(l), H+(l), or Q+(l), respectively, 
for some I £ {0, . . . ,n}, has dimension I 2 , or '(2/ — 1), respectively, and consists of matrices of 

rank not exceeding I. The proof of this assertion is similar for all three cases, for the quaternionic case 
we refer the reader to Proposition IA.6I in the appendix. 

Denote the spaces 5(2), H(2), Q k (2), Q(2) by E 3 ,E 4 ,E 5 ,E 6 , and the cones S+(2), H + (2), Q%{2), 
Q+(2) by A3, A4, A5, Kq, respectively. The index denotes the real dimension of the corresponding 
space or cone. Let T m : E m — > E m , m = 3,4,5,6, be the matrix transposition, or equivalently the 
complex or quaternion conjugation in E m . Note that T m is in Aut(A m ). 

We have the inclusions E m C E n and K m C A„ for m < n. For n £ N+ we have E^^S(n) C S(2n), 
E 4 (g>S(n) C H(2n), E 5 (g>S(n) C Q k {2n) and E 6 (g>S{n) C Q(2n). The space E m ®S(n), m = 3, . . . ,6, 
consists of matrices composed of 4 symmetric n x n blocks. 

Definition 2.7. Let m,n £ N+ and let Abe a mn x ran matrix, with real, complex or quaternionic 
entries. Partition iinraxm blocks A a p (a, /? = !,..., m) of size n x n. Then the partial transpose 
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of A, denoted by A T , will be denned as the result of exchanging the off-diagonal blocks A a p and Ap 
for all a ^ (3. 



A = 



( An A 12 

A21 A22 



\ A m i A m 2 



A lm \ 

A2m 



Amm J 



A = 



( An A 2 i 

A\2 A22 



lni ^2m 



A m 2 



Amm J 



Note that if A is hermitian, i.e. A = A* , then so is A r . 



Definition 2.8. Let A be a hermitian mn x mn matrix, with real, complex or quaternionic entries. 
Then A is said to fulfill the PPT condition or to be a PPT matrix if both A and A T are positive 
semidefinite. 

The set of PPT matrices in the spaces S(m) ($ S(n), Tt(m) (£> S(n), Q(m) <£> S(n), Ti.(m) <g> Tt(n) for 
fixed m,n S N+ is a regular convex cone. 

It is well-known that for matrices in TL{mn) being PPT is a necessary condition for H + (m)<gi H+(n)- 
separability [5]. We can generalize this result in the following way. 

Proposition 2.9. Let V m be one of the matrix spaces S(m) 1 Ti(m), Q(m) and V n one of the spaces 
S(n),TH.(n), Q(n), such that the elements of V m and V n commute. Let K + {m),K + (n) be the corre- 
sponding positive matrix cones. Then any K + (m) (g> K + (n)- separable matrix in V m ® V n is PSD. 

Assume further that K+(m) is invariant with respect to transposition. Then any K + (m) ® K + (n)- 
separable matrix in V rn ®V n is a PPT matrix. 

Proof. It is sufficient to prove the assertions for the extreme rays of the K + {m) (g> -ftf+(n)-separable 
cone. Let the matrix A generate such an extreme ray. Then A can be written as Kronecker product 
A m ® A n , where A m = xx* S K + {m), A n = yy* £ K + (n) generate extreme rays of the corresponding 
PSD cones, and x, y are appropriate column vectors of size to, n, respectively. Hence A = (x(3y)(x<g>y)* 
is PSD, which proves the first part of the proposition. 

Let us prove the second part. We have A T — A m eg) A n . By assumption of the proposition A m 
also generates an extreme ray of K + (m) and can hence be expressed as A m = x'x'* , where x' is an 
appropriate column vector of length m. Then A r = (x' <g> y){x' ® y)* is also PSD and A is a PPT 
matrix. □ 

In particular, the PPT property is necessary for S+(m) <g> S + (n)-, H + (m) <g> S + (n)-, and H + (m) ® 
H + (n)-separability for any to, rt G N+, and for Q + (m) (8> S+(n)-separability for m < 2 and n > 1 
arbitrary (cf. Corollary I A . 8 1 in the appendix). 

Note also that if V n — S(n), then the operation of partial transposition is equivalent to full trans- 
position. If in addition, transposition preserves positivity of matrices in V m <g> V n , then a matrix has 
the PPT property if and only if it is PSD. In particular, this holds for the spaces S(m) ® S(n) and 
H{m)®S{n). 

Let now 3 < to < 6 and n > 1. Denote by S mj „ the cone of K m <B> S+ (n)-separable matrices and 
by T m .n the cone of PPT matrices in the space E m ® S(n). Observe that K m = E m ,i = r. m .i for all 
to. Proposition 12.91 yields the following result. 

□ 

it is a subgroup 



Corollary 2.10. For any 3 < to < 6 and n G N+ 7 we have the inclusion T, m , n C T m ^ n 
Let Qm.n be the group \g®g' \ g € Aut(K m ), g' € Aut(5+(n))}. By Proposition 



of Aut(E TOjn ). Note that the operator T„ 
complex or quaternion conjugation. 

Let us define isomorphisms T m : M. m - 



® id5(„) of partial transposition is in Q„ 

Em, 3 < TO < 6. 



and amounts to 



£ + xi x 2 
X2 x — X\ 

Xq + X\ X2 + IX3 

x 2 - ix 3 x - xi 

Xq + Xi X2+ iX3 + JXA 

X2 — iX3 — jX4, Xo — Xx 

x + xi X2 + ix 3 + jx4 + kx 5 

X2 — IX3 — jx 4 — kx 5 Xq — Xl 



1 3 : {x ,xi,X2) T >-> 
I 4 : (x 0) x 1 ,x 2 ,x 3 ) T h-> 

1 5 : {xo,xi,X2,x 3 ,X4) T i-> 

1 6 : {x ,Xi,X2,X 3 ,X4,X 5 ) T 

It is not hard to check the following result (cf. Corollary IA. 81 in the appendix) 



(1) 
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Lemma 2.11. K m = T m [L m ] for all m = 3, 4, 5, 6. □ 

As a consequence, the group Aut(K m ) is isomorphic to the automorphism group of the Lorentz 
cone L m . 

By virtue of (fTJ) the map I m ^ids( n ) is an isomorphism between the spaces M. m ®S(n) and E m ®S(n). 
Hence we can represent any element of E m ® S (n) in a unique way as an image (T m ® ids(„) ) (X^o 1 e ' ® 
where B , , . . , B OT _i e 5(n). 

Definition 2.12. For any element B = (T m ® id5(„))(^^g 1 e; ® -B;) £ I? m (8 <S(ti), the matrices 
Bq, . . . , B m ^i £ S(n) will be called the components of B. 

Proposition 2.13. The cone r mj „ is invariant with respect to the action of Q m>n , i.e. Q m .n C 

Aut(r mj „). 

Before proceeding to the proof of this proposition, we define j m , m = 3, . . . , 6 to be the set of 
all appropriate matrices S such that the mapping A t— > SAS* is an automorphism of the space E m . 
More precisely, define 73 = GL2 (1^0,74 = GL2(C),jq — GL,2(M) and let 75 be the set of all matrices 
5 £ GL2(H) such that SAS* £ E$ whenever i e £5. The set 75 is a matrix group, its Lie algebra 
given by all quaternionic 2x2 matrices which: i) have a trace with zero i- and j-components, ii) the 
^-components of the off-diagonal elements are zero, and iii) the ^-components of the diagonal elements 
are equal. 

Observe that mappings of the form A 1— » SAS* preserve the PSD matrix cone K m . Let then 
H? n : 7 m — > Aut(if m ) be the group homomorphism assigning to any matrix S £ 7 m the automorphism 
H? n (S) ■ A 1— > SAS* . Denote the image of by G m . Hence G m is the subgroup of all automorphisms 
of E m of the form A i— > SAS* . Let also : GL n (R) — > Aut(5+(n)) be the group homomorphism as- 
signing to any matrix S £ GL n (M) the automorphism H^(S) : A i— > SAS T . This group homomorphism 
is surjective. 

Proof of Proposition \2.1Sl Assume the notations of the proposition. 

The group Aut(L m ) is a Lie group of dimension m ( T ^~ 1 ) _|_ \ anc [ consists of two connected compo- 
nents, corresponding to automorphisms with positive and negative determinants, respectively. Hence 
Aut(A m ) is generated by the connected component of its neutral element and a single automorphism 
with negative determinant. By computing the rank of the Lie algebra of G m one easily determines that 
the dimension of G m equals that of Aut(L m ). It follows [12] that G m contains the connected component 
of the neutral element of A\xt{K m ). For m = 3, 5 the group G m contains automorphisms with negative 
determinants and hence equals Aut(K m ). For m = 4, 6 G m is connected, because 7 m is connected. 
However, in this case the matrix transposition T m has negative determinant as an automorphism of 
E m . Thus Aut(K m ) is generated by G m and T m for all m = 3, . . . , 6. 

Therefore Q m ,n is generated by the following elements. First, elements of the form H^iS) (g> ids(„), 
where S £ 7 m ; second, the element T m <g) ids(„); and third, elements of the form id_B m ®H^{S) 7 where 
S € GL n (R). Let us now consider the action of these generators on the cone L TOjn . Let K + C E m ®S(n) 
be the cone of PSD matrices in E m ® S(n). 

Let S £ 7 m . Then T^ 1 o H? n (S) o T m is an element of G m , and there exists a matrix S' £ ~/ m 
such that T m x o H^{S) o T m = HJ n {S'). The element Hg(S) ® id 5(n) acts on E m <g> 5(n) asin 
(5'(8)/„) J 4(S'(8'/ n )* and hence preserves the cone i4T+. Moreover, we have (iJ^(S')®id5(„))o(T rn ®ids n ) = 
(fq;(S) o T m ) <g> id 5 „ = (T m o H^(S')) ® id 5 „ = (T m ® id s „) o <g> id s(n) ). If we denote the set 

{A r \A£ K+} by Kl, then we get (i/ m (5) <g> id s(n) )[^] = (Hg(S) <8> id s(n) ) o (T m ® ids B )[Jf + ] = 
(T m <g> id s J o Cff^S') (8 id s(n) )[if + ] = (T m ® id 5 „)[A+] = Hence H&(S) (8 id s(n) preserves also 
the cone and therefore also r TO: „ = if + n K+. 

The automorphism (T m (8 id^) is the operator of partial transposition on E m <g> ids( n ) and hence 
preserves T m ^ n by definition. 

Let now S £ GL n (R). Then id Em ®H%(S) acts on E m <g> like A ^ (7 2 ® S)A(7 2 O 5) T 
and hence preserves K + . Moreover, id_E m ®H^{S) commutes with the operator T m ® ids( n ) of partial 
transposition and hence preserves also K+. Therefore it preserves r m>n . 

Thus all generators of G m ,n preserve the cone r min , and so do all other elements. □ 

Proposition 12.131 allows us, when looking for elements in r m rl \ S mj „, to restrict the consideration 
to elements that are in some canonical form with respect to the action of the symmetry group Q m .n, 
since the inclusions in S m ,„ or T m ^ n hold or do not hold for all elements of an orbit simultaneously. 

3 Relations between different cones and trivial cases 

The next three sections aim at proving that the cones Sg^ and are equal. 
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Proposition 3.1. Let 3 < m! < m < 6 and n > n! > and suppose that r mi „ = X m ,„. Then 
V , , — Y. , , 

Proof. Assume the conditions of the proposition. It is sufficient to prove the assertion for n' = n — 
1, m! = m and n' = n, m! = m — 1. 

Let first n' = n — 1, ml = m. Define the linear mapping L„/ „ : S(n') — > S(n) as follows. For 
S G 5(n'), let L n i n (S) be the matrix which has its upper left n' x n' submatrix equal to S and whose 
all other elements are zero. It is not hard to see that if a matrix A £ E m <g S(n) is in the image 
of ids m ®L n :^ n , then A £ S m ,„ if and only if A £ (ids m <g£n',ri)pm,n'] and A £ r m „ if and only if 
A £ (id Bm <g L n ' ,n)[^m,n']- Injectivity of L n ^ n now yields the desired result. 

Let now n = n, m = m — 1 and let C £ T m > Then we have also C £ F mj „ and by assumption 
C £ Tj m>n , Therefore we can represent C as a sum J2b=i M ® Bi, where Ai £ K m and Bi £ S+(n). 
Let now 7r m : R m — + M m be the projection that assigns to any vector (xq, . . . , a; m _2, x m _i) T € R m the 
vector (xo, . . . , x m -2) T £ R m and let 7r T ^ = I m i o-K m oT~^ : E m — > E 1 ,,,,'. Since 7r m [L m ] = we have 
ir^Km] = K m i . Moreover, the restriction of 7r^ to S m / is just idfi m , and the restriction of 7r^ (g <S(n) 
to E m , <g 5(n) is just id Em , 0S(n) . Therefore C = (tt^ <g 5(n))(C) = Ejli ^ ® A with ^ = 
But is an element of K m ', because Ai £ K m . Therefore C is K m > (g S+ (n)-separable and hence in 

^m',n- f— 1 

As one should expect, the equality r„ vn = E m>n is thus easier to prove for smaller n, m. 
By a similar reasoning we can prove the following results. 

Lemma 3.2. Let n,m > 3. TTien t/ie positivity of a matrix in H(m) (g <5(n) is not sufficient for 
H+(m) <g) S+(ri)- separability. The PPT condition in Q(m)<S)S(n) is not sufficient for Q+{m) (g S+(n)- 
separability. 

Proof. Let n, m > 3. Then the cone of positive semidefinite matrices in the space S(m) (g S(n) does 
not coincide with the cone of S+(m) (g S'+(n)-separable matrices. This is a consequence of the fact 
that not every nonnegative definite biquadratic form F(x,y), where x £ M. m , y £ W 1 are two vectors 
of variables, is representable as a sum of squares of bilinear forms [10]. Let then C £ S{m) ® S(n) be 
a PSD matrix which is not 5+(m) <g 5+(n)-separable. Then C is also PSD if considered as an element 
of Ti(m) (g <S(n) and fulfills the PPT condition if considered as an element of Q(m) (g <S(n). 

Now suppose that C is i?+(m)(g)5+(n)-separable. Then there exists an integer N £ N and matrices 
Ai, . . . , Ajy e H + (m), B\, ... , B^ £ S+(n) such that C = M ® Since C is real, we then 

have C = fie(^) ® B t . But i?e(y4;) G ^(m) for all I. Therefore C is S+(m) <g> S' + (n)-separable, 

contrary to our assumption. 

In the same way one shows that Q+ (m) ® S+ (n)-separability of C leads to a contradiction. □ 

It now happens that an element of T, m ^ n can be reduced to elements of cones in lower dimensions. 
We formalize this in the following definitions. 

Definition 3.3. Let 3 < m < 6 and n > 1. Let B £ E m <E)S(n) have components Bq, . . . , B m _i £ S(n). 
If there exists a matrix S £ GL n (M) and positive integers ni, ni with n\+U2 = n such that the matrices 
SBiS T are block-diagonal with blocks B}, Bf of sizes n\ x ni,n<x x ri2 for all / = 0, . . . , m — 1, then 
we call the element B decomposable. We call the elements B 1 = (T m (g 5(ni))(^^ =0 1 e; (g i?; 1 ), 

B 2 = (T m (g 5(n2))(^;ILo 1 e ' ® ^f) ^ oc ^ components of B. 

An element of E m <g 5(n) may have several decompositions and its block components are not 
uniquely defined. 

Definition 3.4. Let 4 < m < 6 and n > 1. If an element B G £J m (g iS(n) has linearly dependent 
components, then we call it reducible. If there exists A £ Aut(J<" m ) such that B' = (A (g id l s(„))(B) G 
E m -i (g S(n), then we call B' a reduction of B. 

Lemma 3.5. Lei B G r m n be reducible. Then B /ias a reduction. 

Proof. Assume the conditions of the lemma and let Bq, . . . , S m -i be the components of B. Then there 
exists a nonzero vector A = (Ao, . . . , A m _i) T G K m such that X^o 1 ^i^i = 0. 

If B = 0, then B G £? m _i(g<S(n) and B is a reduction of itself. Hence let us assume that B ^ 0. Since 
B £ r m „, the diagonal elements of the components Bq, . . . , -B m -i cannot all be zero simultaneously. 
Without restriction of generality let us assume that the (1, l)-elements of Bq, ... , £? m _i form a nonzero 
vector b = {b , . . . ,6 m _i) T £ R m . Define b' = (h, . . . , 6 m _i) T , A' = (A 1; . . . , X m -i) T G E m_1 . 

The element b = T m (b) £ E m is a principal 2x2 submatrix of B. Since B G T mjn , we have b y 
and hence b £ K m , b £ L m , and ^ 6o > ||o'||. Moreover, (6, A) = 6oAo + (b',X') = 0. It follows that 
&o | Aq | < 1 1 6' 1 1 1 1 A' 1 1 < &o||A'|| and |Aq| < ||A'||. We distinguish two cases. 
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i) | Ao | = 1 1 A' 1 1 7^ 0. Let without restriction of generality Ao = 1. Now choose an orthogonal 
(to— 1) x (to— 1) matrix U such that its first row is given by A' . Define the automorphism A' £ Aut(L m ) 
by A' — diag(l,J7) and let A = X m A'T^- be the corresponding automorphism of K m . Denote the 
components of B' = (A (8 id lS ( n - ) )(B) by B' , . . . , B' m _ 1 . We now compute B' and B[. We have 
B' = Bo, B[ = jyi^i 1 A;S; = — Ao-Bo- Hence B' + B[ = 0. Note that this sum is the upper left n x n 
block of the matrix B'. But A® ids(„) £ Gm,n, hence B' £ T m ^ n by Proposition 12.131 and in particular 
B' is PSD. It follows that the first n rows and columns of B' are zero. In particular, we get B' m _ x = 0, 
and B' £ E m _ x ® S(n). 

ii) | Ao | < ||A'||. Let without restriction of generality ||A'|| 2 = Aq + 1. Then there exists a number 
£ and a unit length vector v £ R m_1 such that Ao = sinh£, \' T — cosh£u T . Let now A\ £ Aut(Z m ) 
be a hyperbolic rotation in the (eo, e m _i)-plane by the angle £, and let A' 2 £ Aut(L m ) be a rotation 
in the linear subspace spanned by {ei, . . . , e TO _i}, given by an orthogonal (m — 1) x (m — 1) matrix 
which has v as its last row. The last row of the matrix A' = A\A^ <= Aut(i m ) is then given by 
(sinh£, cosh£ V\,..., cosh£u m -i) = A T . Define now the automorphism A = I m A'I^ £ Ant(K m ) and 
let B' = (A ® id,5(„))(B). The last component of B' is then given by B' m _ 1 = YmL^ M^l = 0- Hence 
B' £ E m -i ® S(n). □ 

We can hence reduce decomposable elements to elements in spaces with smaller n, and reducible 
elements to elements in spaces with smaller to. By using Proposition 12.131 and applying the lines of 
reasoning in the proof of Proposition 13.11 we arrive at the following proposition. 

Proposition 3.6. A decomposable element ofT m n is in S TO „ if and only if its block components are 
in E m iu and S m ll2 , respectively. A reducible element ofT mn is in ^> m ,n if and only if its reduction 
is in E m _ li7l . □ 

Corollary 3.7. Suppose that r m ,„ = S m)W for some n,m. Then all decomposable elements ofY m n+ i 
are in S m>n -|_i and all reducible elements o/Tm+Ln are in S m +i, n . 

Let us now consider the cases to = 3 and n = 2. 
Theorem 3.8. „ = £3^ for any n > 1. 

Proof. The cone r3 jra is the cone of real symmetric PSD block-Hankel matrices of size 2n x 2n. However, 
such matrices are known to be separable. This follows from the spectral factorization theorem for 
quadratic matrix- valued polynomials in one variable [16] . □ 

By a similar reasoning applied to complex hermitian block-Hankel matrices we obtain the following 
result. 

Theorem 3.9. Let n > 1. A matrix in 5(2) ® TL(n) is PSD if and only if it is S + (2) eg) H + (n)- 
separable. □ 

Theorem 3.10. r m ,2 = S m ,2 for to = 3, . . . , 6. 

Proof. The space 5(2) is 3-dimensional. Hence any quadruple of real symmetric 2x2 matrices is 
linearly dependent, and for to > 4 the space E m ® S(2) consists of reducible elements only. Since 
T3,2 = £3,2 by the previous theorem, we have r m ,2 = £m,2 for arbitrary to = 3, . . . , 6 by repeated 
application of Corollary 13. 71 □ 

Lemma 3.11. Let B £ be partitioned in symmetric 3x3 blocks as follows 

B = ( — Bl2 \ (2) 

\ B\2 B 2 2 J 

Then either B is reducible or decomposable, and hence in £4,3; or in the orbit o/B with respect to the 
action of the group ^4,3 there exists an element B' such that its upper left 3x3 subblock B[ x equals 
h- 

Proof. Let the assumptions of the lemma hold. Then we have B 0. Hence if a real vector v £ M 3 is 
in the kernel of one of the matrices Bn, B22, then it is also in the kernels of B12, B\2- 
Suppose that neither B\\ nor -B22 are of full rank. 

If the kernels of B\\,B22 have a nontrivial intersection, then B is decomposable, because this 
intersection will be in the kernel of all four components of B. 

If the intersection of the kernels is trivial, then the kernel of B12 contains two linearly independent 
real vectors. Then the real and imaginary parts of -B12 must be linearly dependent, because they are 
symmetric, sharing a 2-dimensional kernel and of rank not exceeding 1. Hence B is reducible. 
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By Corollary 13.71 and Theorems 13.81 and 13.101 any reducible or decomposable element of is in 
£4,3- 

Let us now suppose that at least one of the matrices B\ 1,-622 has full rank and is hence positive 
definite (PD). We can assume without loss of generality that this is Bn, otherwise we pass to a matrix 
in the same orbit by first applying the automorphism ifj (cr2)<S>ids(3) G ^4,3, where 02 is the non-trivial 
2x2 permutation matrix. Then the matrix B' = (ide 4 (g)i?^((i3ii) _1 / 2 ))(B) is in the orbit of B and 
we have B' xl = I 3 . □ 

Theorem 3.12. r 4 , 3 = £4,3. 

Proof. We shall show that any extreme ray of is in £4,3. 

Let B generate an extreme ray of T4.3. Note that T4.„ is isomorphic to the intersection of the 
4n 2 -dimensional cone H + (2n) with the 2n(n + l)-dimensional subspace of matrices consisting of four 
symmetric n x n blocks. Then the face of B in H+(6) has at most dimension 36 — 24+ 1 = 13 by 
Lemma \2. 61 Hence the PSD matrix B G H+(6) has at most rank 3. 

By Lemma 13.111 we can assume without restriction of generality that B is partitioned as in |(3J) with 
Bn = I 3 . Since B has rank 3, we can factorize it as 

( h_ B 12 \ ( I 3 \( h V 
V B 12 B 22 ) \W )\W ) 

with W = lh 2 = W T . It follows that WW = WW* = B 22 = (WW*) T = WW. Therefore the 
real and imaginary parts of W are symmetric and commute. In particular they can be simultaneously 
diagonalized by applying an orthogonal transformation U. Applying the automorphism id_g 4 ®H^(U) 
to B, we thus simultaneously diagonalize all four matrices Bn, B\ 2 , B\ 2 , B 22 , and B is decomposable 
and hence separable. 

We have proven that all extreme rays of T4 j 3 are in £4,3. But then we have r4,3 = £4,3. □ 

In this section we have investigated the relationship between the cones E m , n ,r mi „ for different 
dimensions to, n. We have defined two properties of elements in E m <X> S(n), namely those of being 
decomposable and reducible. Our next goal is to show that r mj „ C £ m ,n for n = 3, m = 5, 6 (the 
converse inclusion being trivial). We have shown in this section that this inclusion is valid for elements 
possessing the above-cited properties, provided the relation T m>n — £ m , n holds for the respective cones 
of smaller dimension. We have proven this relation for n < 2; m < 3; and m = 4, n = 3. This allows 
us to concentrate on non-decomposable and non-reducible elements in the proofs of the main results 
in the next two sections, which essentially amounts to imposing certain non-degeneracy conditions. 



4 r 5;3 — E 5;3 

The structure of the proof of this equality resembles that of the proof of Theorem 13. 121 We show that 
every element of ^3 generating an extreme ray is in £5,3. 

First we derive some properties of real 3x3 matrices. Let A(n) be the space of real skew-symmetric 
matrices of size n x n. The space .4.(3) is isomorphic to M 3 . We define an isomorphism V : -4(3) — > R 3 
by 

(0 a i2 -a 3 i \ / a 23 \ 
-012 a 23 J !— ► I 0,31 J . 
031 ~«23 / \ a 12 J 

If v, w G R 3 are column vectors, then this isomorphism maps the skew-symmetric matrix vw T — wv T 
to the cross-product v X w. 

Define now a group homomorphism H Q : GL 3 (R) GL 3 (R) by [H®(C)](v) = V(C[V- 1 («)]C T ) 
for any v e R 3 , C € GL 3 (M). Note that iJ Q (-C*) = iJ Q (C) for all C € GL 3 (R). One also easily 
checks that the induced Lie algebra homomorphism has a trivial kernel. Therefore the image of is 
the connected component of the identity matrix ^3 and consists of the real 3x3 matrices with positive 
determinant. Direct calculus shows that H®(C) = (det C)C~ T . 

For 3 column vectors u,v,w G R 3 , let (u,v,w) be the 3x3 matrix composed of these column 
vectors. The following result can be checked by direct computation. 

Lemma 4.1. For any three vectors u,v,w G R 3 we have (u x v) x (v x w) — v det(u, v, w). Moreover, 
the following assertions are equivalent. 

i) det(u, v, w) = 0, 

ii) u,v,w are linearly dependent, 



S 



Hi) u x v, v x w, w x u are linearly dependent, 
iv)uxv,vxw,wxu are all proportional. 



□ 



Corollary 4.2. Let u,v,w € R 3 be linearly independent vectors. Then dct(v x u, w x v, u x w) < 0. 
Proof. By repeated application of the previous lemma we get 

v x w = [(it x t>) x (i> x w) dct _1 (u, v, w)] x [(v x to) x x u) det~ 1 (u, w, u)] 
= det _2 (w, v, w)[(u x v) x (v x w)] x [(w x u) x (id x u)] 
= det _2 (w, v, w)[(v x w) det(u x v, v x w, w x u)]. 

It follows that det(w x it, w x u, it x w) = — det(u xv,vxw,wxu) = — det 2 (it, v, w) < 0. □ 

We now investigate 2-dimensional traceless linear subspaces C in 5(3). To each such subspace, 
we will assign a sign a(£) € { — 1,0, +1} in the following way. Let {Si, S2, S3, I3} be a basis of 
the orthogonal complement £-*- of £ Consider the vectors V12 = V(SiS2 — S2Si),V23 = V(S , 2<S I 3 — 
SsS2),V3\ — V(SsSi — Si S3) G M 3 . Now define the sign of C as a(£) = sgndet(w23, W31, ^12). 

Lemma 4.3. The sign o~(L) is well-defined, i.e. does not depend on the choice of the basis of C 1 - . 

Proof. Let us first remark that C 1 - has dimension 4 and contains the identity matrix I 3 , hence the 
claimed choice of its basis is possible. 

Let now {Si, S2, S3, 13} and {S[, S' 2 , S' 3 , 13} be two bases of C 1 - . Observe that [A + aI n , B + f3I n ] = 
[A, B] for any two nxn matrices and any real scalars a, (3. Since <j{£) depends only on the pairwise 
commutators of the basis elements Si,S[, we can assume without loss of generality that these matrices 
are traceless. Note that then both Si and S[ span the same 3-dimcnsional space, namely the orthogonal 
complement of C in the subspace of traceless symmetric matrices. We hence find a regular 3x3 
matrix C such that S' t = J2a=i Ci a S a for I = 1,2,3. Here the indexation of C denotes its elements. 
Define v a fj = V([S a , Sp]), v' a/3 = V([S' a , S^]), a,f3 = 1,2,3. Then we have by the bilinearity of the 
matrix commutator that v' a g = YMm=i^ctl v lm^mB f° r a i P — 1 , 2, 3. It follows that (^23,^31,^12) = 
{v23,vzi,vi 2 ){H Q {C)) T = (detC) (w23,«3i,wi2)C~ 1 and det(u 23 , v' 31 , v[ 2 ) = (det C) 2 det(w 2 3, V31, V12). 
Thus the signs of det (^35 ^31 , v' 12 ) an d det(v23, V31, V12) are equal. □ 

Definition 4.4. We call two 2-dimensional traceless subspaces £, CJ C 5(3) equivalent if there exists 
an orthogonal 3x3 matrix U such that £ = {USU T \ S G £}. 

Lemma 4.5. For equivalent subspaces £,£ we have a(£) = <r(£). 

Proof. Let U be the orthogonal matrix realizing the equivalence. Let {Si, S2, S3, 13} be a basis 
of C 1 - and define S[ = USiU T , I = 1,2,3. Then {S[ , S 2 , S 3 , 7 3 } is a basis of £^ . Let further 
v a p = V([S a ,S p ]), v' af} = V([S' a ,S'p}), a, /3 = 1,2,3. Since we have [S' a ,S' ] = U[S a ,S p ]U T for 
all a, [3 = 1,2,3, it follows that v' a0 = [HQ(U)]v a(} and (^3^31^12) = [H Q (U)](v 23 , V31, v 12 ) = 
(det U) U{v 2 3, w 3 i, v l2 ). This finally yields a{£) = a{£). □ 

Let us now consider symmetric tensors S Q( a 7 of order 3 in ]R 3 . Since there are 10 independent 
components, these tensors form a 10-dimensional real vector space. 

Definition 4.6. We shall say that a symmetric tensor S Q , ( a 7 of order 3 satisfies the 5-condition if 
X^=i(Sa/3KS 7T;K - S 7l 3 K S ar]K ) = S 7/3 S ar) - S a 0S ln for all a,/3,7,?7 = 1,2,3, where S is the Kronecker 
symbol (S a p = 1 if a = [3 and 8 a p = otherwise). 

Definition 4.7. Let S Q/ 3 7 be a symmetric tensor. Let the matrix components of S Q /3 7 be three matrices 
S 1 , S 2 , S 3 e 5(3) defined elementwise by S l afj = S afM , a, (3,1 = 1, 2, 3. 

Remark 4.8. The (5-condition is equivalent to the condition 5Zf =1 (det S')(S') _1 = —I3, where S l are 
the matrix components of the tensor. Here the function A 1— » (det A)A _1 is understood to be extended 
by continuity to singular matrices A E 5(3). 

Lemma 4.9. The 5-condition is rotationally invariant, i.e. ifU is an orthogonal 3x3 matrix, then the 
tensor S^ l(37 = Y% <t> £=1 Ua-qUp^U^^S^ satisfies the 5-condition if and only if the tensor S a pj does 
so. 

The proof is by direct calculation using the relation y] ^ Ua^Upj = 8 a p. 

Lemma 4.10. Let S a /3j be a symmetric tensor satisfying the 5-condition with matrix components S l , 
I = 1,2, 3. Then the matrices {S 1 , S 2 , S 3 , 13} are linearly independent. 
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Proof. We proof the lemma from the contrary. Assume the conditions of the lemma and suppose that 
there exist co,ci,C2,C3 £ R, not all equal zero, such that C0/3 + Jw=i c lS l = 0. By Lemma 14.91 we 
can assume without loss of generality that C2 = C3 = 0. Then there exists c ^ such that S* 1 = c/3, 
or Sip-y — cSpj for all /3, 7 = 1,2,3. Since S a p-y satisfies the (5-condition, we have in particular 
^^ =1 (S , h k S , 22k — 5 , 2i k S'i2k) = ^2i<5i2 — <5n^22- But the left-hand side of this equation simplifies to 0, 
whereas the right-hand side simplifies to —1, which leads to a contradiction. □ 

We now come to a result linking the sign of a traceless 2-dimensional subspace £ C 5(3) to 
symmetric tensors satisfying the <5-condition. 

Theorem 4.11. 2-dimensional traceless subspaces of 5(3) with sign —1 are in correspondence with 
symmetric tensors satisfying the 5-condition. Namely, if £ is such a subspace, then there exists a 
symmetric tensor 5 Q /3 7 with matrix components S l , I = 1,2,3 satisfying the 5-condition such that 
{S 1 , S 2 , S 3 , 13} is a basis of the space £ . Conversely, if S a p^ is such a tensor, then the orthogonal 
complement of the linear span of the set {S , S 2 , S 3 , I3} has sign —1. 

Proof. Let e\, e 2 , e^ be the canonical orthonormal basis vectors of R 3 . 

Let S a (3-y be a tensor with matrix components S l , I = 1,2,3 satisfying the (5-condition. Then we 
have S 1 S 2 - S 2 S 1 = e 2 e\ - ei e% and V([S X ,S 2 \) = e 2 x ei = -e 3 . Similarly, V([S 2 ,S 3 }) = -e x and 
V([S 3 , S 1 }) = — e 2 . Let £ be the orthogonal complement of the linear span of the set {S 1 , S 2 , S 3 ,l3}. 
Then by definition o~(£) — sgndet(— e\, — e 2 , —e^) = —1, which proves the second part of the theorem. 

Let us prove the first part. Let £ C 5(3) be a 2-dimensional traceless subspace with sign o~{C) = — 1. 
Any point in 5(3) can be viewed as a homogeneous quadratic form on R 3 , or equivalently, as a quadratic 
map from RP 2 to R. We are interested in the number of points in RP 2 which are mapped to zero by 
all elements of a C. Denote N c = {x e R 3 | x T Sx = V S E £}. 

The determinant as a scalar function on C is odd and hence possesses zeros on £ \ {0}. Since the 
matrices in C are traceless, we can find a matrix in C with eigenvalues —1, 0, +1. By conjugation with 
an appropriate orthogonal matrix U we can transform it to the matrix 

10 
Ni = I 1 


Let £' — {USU T I S £ £}. Then the subspaces £ and £' are equivalent and we have N\ £ £'. A 
generator the orthogonal complement of Ni in £' will be of the form 

for some numbers a,b,c,d S R, which do not equal zero simultaneously. Note that by Lemma \4. 5 1 we 
have <j{£') = — 1. Let us first treat several degenerate cases. 

1. a = b — 0. In this case a basis of £ ,A ~ is given by {S\, S2, S3, 13} with 

/ 1 \ / \ 

Si = , s 2 = I 1 I , s 3 = 

\ / 

We get a{£') = 0, because [Si, S 2 ] = 0. 

2. a = — b ^ 0. In this case a basis of £ ,A ~ is given by {Si, S 2 , S3, 13} with 

-ad -cd \ / -ac a 2 + d 2 

Si = I I , S 2 = I ad a 2 + c 2 , S3 = ac -cd 

/ \ a 2 + d 2 -cd 

We get det(v 23 ,v 31 ,v 12 ) = a 4 (a 2 + c 2 + d 2 ) 2 > and a(£') = 1. Here v a/3 = V([S a ,Sp]). 

Hence these two cases do not satisfy the conditions of the theorem and we can assume a + b ^ 0. 
The set Nc> is given by those vectors x = (xi, x 2 , X3) T £ R 3 that satisfy x T N\X = x T N 2 x = 0. In 
particular, x must satisfy x±x 2 — 0. If we define 

N 21 =( a C , V N 22 = 1 " 















• 





• 














c —a — b ) ' V d —a — b 
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then we get 




No = {\ | | (i x x 3 )iV 21 ( -j=0|u|^ 2 j|(, 2 x 3 )N 22 )=0}. (3) 

Note that these two sets have a trivial intersection due to the condition a + b ^ 0. Let D x = det N 2 i, 
D 2 = det N 22 . We have Di+ D 2 = — (a + b) 2 — c 2 — d 2 < and at least one of these two determinants 
is negative. Let this be D\ without restriction of generality (otherwise we exchange a, 6; c,d; X\,x 2 ; 
and N 2 i,N 22 ). 

A basis of C lA ~ is then given by {S\, S 2 , S 3 , 1%} with 



(a + b) 2 (a + b)(c + ^Dl) 

s t i 

(a + b)(c + ^f^D~i) (c + T^A") 2 

(a + b) 2 (a + 6)(c-y^DT) 







,2 



(a + b)(c- V=DI) (c-y^A" 

—ad — cd 

,V, = i -bd a 2 + ab + b 2 + c 2 

-cd a 2 + ab + b 2 + c 2 (a + b)d 

We then get det («23,f 31,^12) = — 4a6(— D±)(a + b) A (a 2 + ab + b 2 + c 2 ) 2 and (t(C') = - sgnab. Since 
<j{£) = — 1) we get a6 > 0. Since we can multiply N 2 by —1, we can assume without restriction of 
generality that a,b > 0. We have D 2 = —ab — b 2 — d 2 < 0. Hence both N 21 and N 22 have negative 
determinant. But the number of linearly independent vectors in the two subsets on the right-hand side 
of © depends precisely on the sign of these determinants. Namely, if N ^ is a real symmetric 2x2 
matrix with negative determinant, then the equation y T Ny — has two linearly independent solutions 
y E R 2 . Thus Nc contains 4 linear 1-dimensional subspaces of M 3 . 

Let x , I = 1, . . . ,4 be generators of these subspaces, defined by the columns of the matrix 




a-\-b fa a+b 




X = 



Each 3 of these vectors span the whole space R 3 and 

(x a ,x ) = -l V a, (3 = 1,2,3,4, a £ (3. (4) 

The rank 1 matrices x [ (x l ) T <E S(3) are linearly independent and hence form a basis of the space £ ,J ~ . 

Let us define symmetric matrices X, I = 1, .. .,4 by X 1 = x l (x l ) T — I3, and positive numbers 
q = 1+ u.ii2 ■ Consider the symmetric 4x4 matrix X T X. It has rank 3 and is hence singular. Its diagonal 

elements are given by |a;'| 2 , while its off-diagonal elements equal —1 by (j4|). Setting the determinant 
to zero and transforming the obtained relation, we get Ef=i Q = 1. Then using ([4]) one checks that 
(E?=i c & 1 , x p ) = (Eti ciX 1 , X p( X p) t ) = for all p = 1, ... ,4. This implies £? =1 qx' = £? =1 c/X' = 
0. Hence the matrices X 1 span a 3-dimensional subspace C 3 C £ ,J ~, which does not contain the matrix 
J3, and there exists a unique linear mapping X : ]R 3 — > £3 such that A^ar) = X z for all i = 1, . . . , 4. 
Moreover, we have again by © that X"^ = X^x" and X Q X /3 - X$X a = x p {x a ) T - x a {x p ) T for all 
a, 13= 1,2,3,4. 

By bilinearity these relations extend to X(v)w — X(w)v and X(v)X(w) — X{w)X{v) = wv T — vw T 
for all v,w £ K 3 . Let now S l = X(ei), I — 1,2,3. Then by the first relation the matrices S l are the 
matrix components of some symmetric tensor S' af3j , and by the second relation this tensor satisfies the 
(^-condition. Thus S' a p satisfies the assertions of the theorem with respect to the subspace £ . For 

reference, the components of this tensor equal S' in — — 7=, S' 222 = -j=, S' 113 — S' 223 = — 



_ Tab' "222 - "113 - V b' ^223 

^333 = VT — \/a ' au other independent components being zero. 

Finally, define the tensor S a /3-y — J2ri & g=i ^arpj^^j^'^ o w nere U is the orthogonal matrix 
realizing the equivalence of the subspaces £ and £'. By Lemma T4T51 this tensor satisfies the assertions 
of the theorem with respect to £. □ 

Finally we provide the following auxiliary result about quaternion vectors. 
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Lemma 4.12. Let v £ H 3 and p 6 i, p / be given. Then there exists a unit quaternion q such that 
Re(vqp) = 0. Similarly, there exists a unit quaternion q' such that Re(q'vp) = 0. 

Proof. The equation Re(vqp) — (or Re(q'vp) = 0) amounts to 3 linear relations on the 4 real 
components of q (or q'). Hence there exists a non-zero solution, which can be normalized. □ 

In particular, we can multiply any vector in H 3 by a unit quaternion from the left or from the right 
such that a given component (e.g. the real part) vanishes. 

Let us return to the cones and Ss )n . Our goal is to show that these cones are equal for n = 3. 

Lemma 4.13. Any matrix B £ ^5,3 is either reducible or decomposable, and hence in £5,3, or in the 
orbit of B with respect to the action of the group C/5,3 there exists an element B' such that its upper 
left 3x3 block B'n equals I3. □ 
The proof is similar to the proof of Lemma 13.111 with a reference to Theorem 13.121 instead of 
Theorem ESI 



Let now B be the generator of an extreme ray of ^3. Our goal is to show that B £ £5,3. Let us 
apply the dimensional argument Lemma 12.61 



Corollary 4.14. Let B generate an extreme ray ofT^^. Then B has at most rank 4- 

Proof. Any element B £ satisfies B y 0. The space 2(6) has dimension 66, whereas its subspace 
E5 (g> S(n) has dimension 30. Hence the face in Q+(6) of the extreme ray B has at most dimension 
37 by Lemma T2.6I Since the space Q(5) has dimension 45, a rank 5 element cannot be extremal (cf. 
Proposition [A6] in the appendix). □ 

Lemma 4.15. Let B generate an extreme ray ofT^^. Then B £ £5,3. 

Proof. Assume the conditions of the lemma. Let B be partitioned in 4 symmetric n x n blocks as in 
@. By Lemma T4. 131 and Proposition 12. 131 we can assume without loss of generality that B\\ = L3 and 
B is not reducible. 

Denote the 3 x 3 matrix Uj^ by W . Then B 2 2 - WW* is a PSD matrix, which by Corollary [4. 141 
has rank at most 1 (cf. Lemma [A. 71 in the appendix). Hence there exists a quaternionic vector Z £ H 3 
such that ZZ* + WW* = 822- Here Z can also be the zero vector. The vector Z is determined up to 
multiplication with a scalar of unit norm from the right. By Lemma 14.121 we can find a unit quaternion 
q such that Zq has a zero ^-component. Let us hence assume without loss of generality that Z has a 
zero /c-componcnt. 

Let us factorizc B. We have 



B 



h W* \ = I \ I 
^£22/ { w z { w z 



Note that W is symmetric and with zero fc-component. We have as well that WW* + ZZ* = B22 is 
symmetric and real. This implies 

WiW r - W r Wi + ZiZ? - Z r Zf = 0, 

W 3 W r - W r Wj + ZjZ? - z r zj = 0, 

WjWi-WiWj+ZjZT -ZiZf = 0, (5) 

where the indices denote the corresponding components of the matrix W and the vector Z. Let us 
denote the row vector Z T by Y. Then we have Z* = Y. Since B is not reducible, the four matrices 
W r , Wi, Wj, I3 are linearly independent. Denote the orthogonal complement of their linear span by £,. 
We shall consider several cases. 

1. The components of Z are linearly independent. 

Condition and Corollary 14.21 imply that C has negative sign. By Theorem 14.111 there exists 
a symmetric third order tensor S a pj fulfilling the (5-condition whose matrix components span a 3- 
dimensional subspace C3 C . This subspace is transversal to ^3, therefore there exist real numbers 
x r , Xi, Xj such that the matrices W r — x r L 3 , Wi — Xi/3, Wj — xjl 3 span £ 3 . 

Let us now consider the symmetric completion of the 3x4 quaternionic matrix (W Z) to the square 
matrix 

N=(™ Z 
Y x 
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where x = x r + ixi + jxj is a quaternion with zero fc-component. The matrix TV has zero fc-component 
and mutually commuting real, i- and j-component. Namely, the commutation condition amounts to 
the equations 

[W a ,Wp] = ZpY a -Z a Y , 
W a Zf3 + Z a Xfj = WpZ a + Zpx a , 
Y a Zp + X a Xfj = YpZ a + xpx a , 

for every pair of indices a, /3 from {r,i,j}. The first set of equations is precisely (O, the last set of 
equations is satisfied for any x because Y a = Z^ for all indices a. These second set follows from the 
symmetricity of the tensor S a 0-y . 

Since the matrices N r ,Ni,Nj mutually commute and are symmetric, they share a common or- 
thonormal set of real eigenvectors v±, t>2, 1)3, V4 £ K 4 . This set is also a set of eigenvectors to the 
quaternionic matrix N, with eigenvalues 51,52,93,94 that have zero fc-component. Hence we have 
(W Z)vi = qi{I% 0)vi = qivi for all I £ {1,2,3,4}. Here vi £ M 3 consists of the first three components 
of vi. If we decompose the identity matrix I4 as Yli=i v i v I \ we can rewrite B as 

This matrix is thus K 5 (g> S' + (n)-separable, and B £ £5,3- 
2. The components of Z are linearly dependent. 

Let z £ R 3 be a nonzero vector in the orthogonal complement to the subspace spanned by the 
components of Z. Then by (|5|) we have [W r , Wi]z = [Wi, Wj]z = [Wj, W r ]z — 0. By bilinearity we then 
have [U, V]z = for any two matrices U, V 6 £- L . Define the linear subspace Sw = {U E \ Uz = 0} 
and let S z C M 3 be the linear span of the vectors z = I3Z, W r z, WiZ, WjZ. Obviously we have 

dim Sw + dimS z = 4. (6) 

For any matrices U £ Sw and V £ C 1 - we have U(Vz) = V(Uz) = 0. It follows that for any matrix 
U £ Sw and any vector y £ S z we have Uy = 0. If now dim<S 2 = 3, then Sw can consist only of the 
zero matrix, which contradicts ©• Likewise, if dim5j = 2, then all matrices in Sw are proportional 
to the rank 1 matrix vv T , where v is the orthogonal complement to S z . Hence dim Sw < 1, which also 
contradicts ([5]). Therefore dimS z = 1, and z is a common eigenvector of all matrices in C . It follows 
that z is an eigenvector of W,W*, and also ZZ*, because Z* z — 0. Thus z is a common eigenvector 
of the 5 matrices B\\ = Ia,W r ,Wi,Wj, -B22 and these matrices are linearly dependent. But then the 
components of B are also linearly dependent, which contradicts the non-reducibility of B. 

This completes the proof of Lemma 14.151 □ 

We obtain the following theorem. 

Theorem 4.16. F 5j3 = E 5j3 . □ 



5 r6,3 = ^6,3 

In this section we prove the analogue of Theorem 14.161 for dimension m — 6, following essentially the 
same line of reasoning. First we provide some auxiliary results. 

Lemma 5.1. Let v, w £ H 3 be quaternionic column vectors of length 3 such that Revv* = Reww*. 
Then there exist unit norm quaternions h,h' such that h'v = wh. 

Proof. Let the assumptions of the lemma hold. The quaternions H can be considered as a 4-dimensional 
real vector space, equipped with the Euclidean scalar product (a, b) = Re ab. Let W : H — > K 4 be the 
corresponding isomorphism. We have for any pair of indices a, (3 — 1,2,3 that Rev a vJ = Rew a Wp- 
We then get (W(v a ), W{vp)) = (W(iv a ), W(wp)). Therefore there exists an orthogonal 4 x 4-matrix 
U such that W(w a ) — UW(v a ) for all a = 1, 2, 3. Since v, w have only 3 components, this matrix can 
be chosen to have determinant 1. But the group generated by multiplication of H by unit quaternions 
from the left and from the right is known to be isomorphic to the special orthogonal group 50(4). 
This can be easily checked by comparing the Lie algebras of these groups. Hence there exist unit norm 
quaternions h, h! such that w = h'vh" 1 . □ 

We provide the following lemma on the characteristic polynomial of a 3 x 3 hermitian quaternionic 
matrix. 
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Lemma 5.2. Let A be a hermitian quaternionic 3x3 matrix with elements a a p — a 1 a p + ia l a p + 3 ' al a p + 
ka/^g, a, (3 — 1,2,3. Then the coefficient c\ of its characteristic polynomial A 3 + C2A 2 + c\\ + cq is 
given by 

c l = -(4s) 2 + a ll a 22 " (als) 2 + a 22 a 33 + a ll a 33 ~ ( a 12) 2 ~ K 3 ) 2 ~ Ks) 2 ~ ( a 2sf 

-(«i 2 ) 2 - fe) 2 - K2) 2 - (af 3 ) 2 - (<4) 2 - (a* 3 ) 2 . □ 

The proof is by direct calculation (cf. Proposition IA.3I in the appendix) . 

The partial transposition in the space E§ ® S(n) amounts to quaternionic conjugation. Hence the 
cone is composed of those B 6 Eg <g) 5(3) that satisfy the linear matrix inequalities 

Bn B12 \ =■ / Bn B12 



Here B is partitioned in 4 symmetric 3x3 blocks as in ([2|). 

Lemma 5.3. Let B G r 6 ,3. Then either B is reducible or decomposable, and hence in S 6j3 , or in the 
orbit o/B with respect to the action of the group Q 6 3 there exists an element B' whose upper left 3x3 
block B'n equals L^. □ 
The proof is similar to the proof of Lemma 13. Ill with a reference to Theorem 14.161 instead of 
Theorem 13.. 



We now proceed as in the previous section and look for elements B € r 6 3 which are not in S 6 3 . 
Let B be partitioned as in ||7J). By Lemma [5751 and Proposition 12.131 we can assume without loss of 
generality that B is not reducible and that Bn = L3. 

We can factorize the PSD matrices B, B as 

n=f 1 °U J °V n=(L \(L V 

\ W Z J \ W Z J ' \ W Y* J \ W Y* J ' { ' 

where W = Byi , and Y, Z are quaternionic matrices of appropriate size (cf. Lemma IA.7I in the ap- 
pendix). Note that W and WW + ZZ* = WW + Y*Y = B 2 2 are symmetric matrices. 

As in the previous section, denote the components of W, Z, Y by adding a corresponding index. 
The imaginary part of WW + ZZ* is zero, which yields 

[WuW^ + lW^W^+ZiZj -Z r Zj + Z k Zj -Z^Zl = 0, 
[W 3 ,W r ] + [Wi,W k ] + ZjZ? - Z r Zj + ZiZl - z k zj = 0, 
[W k ,W r ] + [W h Wi]+Z k Z^ -Z r Zl + Z j Zj -ZiZj - 0. (9) 

In a similar way, setting the imaginary part of WW + Y*Y to zero, we obtain 

- [Wi , W r ] + [W k ,Wj]- Y?Y r + YjYi + Y k T Yj - YjY k = 0, 
- [Wj , W r ] + [Wi,W k \- YjY r + Y r T Y 3 + Y?Y k - Y k T Y t - 0, 
-[W^W^ + lW^W^-Y^Yr+Y^Yk + Y^Y.-Y^Y, = 0. (10) 



Let us apply the dimensional argument Lemma 12.61 Referring to the notations of this lemma, we 
choose K to be the 132-dimensional direct product Q+(6) x Q + (6). Its faces are isomorphic to direct 
products Q + (li) x Q + (l 2 ) with l\ < 6, l 2 < 6 (cf. Proposition IA.6I in the appendix). The subspace L 
consists of all pairs (B, B) E <2(6) x Q(6) with B consisting of four symmetric 3x3 blocks. This space 
is 36-dimcnsional. The intersection K' = K n L is isomorphic to the cone of blockwise symmetric PPT 
matrices in Q + (6), i.e. to r6,3- By Lemma |2.6[ an extremal ray of r6,3 lies in a face of K that has at 
most dimension 97. Now note that the space <2(5) x Q(6) has dimension 45 + 66 = 111. Hence if B or 
B are of full rank, then the corresponding conjugate can be at most of rank 4. We get the following 
corollary. 

Corollary 5.4. Let B, factorized as in (B\), define an extreme ray ofT^^. Then either rk ZZ* < 1 or 
rk Y*Y < 1 or rk ZZ* = rk Y*Y = 2. □ 



We treat these cases separately. Let the assumptions of Corollary 15.41 hold and suppose that B is 
not reducible. 

1. The matrix ZZ* has at most rank 1, i.e. Z is a quaternionic vector, which can be zero. 
1.1. The four components of Z lie in a 2-dimensional linear subspace of R 3 . 
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Then there exists a non-zero vector z G R 3 such that Z* z = 0. We have WW — WW = 
W(2W r — W) — (2W r — W)W = 2[W, W r ]. It follows that ||Fz|| 2 = z*Y*Y z = z* {2[W,W r ] + Z Z*)z = 
2z T [W,W r ]z = 0. The last equality follows from the fact that z is real and hence commutes with 
quaternions, and that [W, W r ] is skew-symmetric. Therefore we also have Yz = 0. Multiplying © and 
pop by z from the right and combining, we obtain [W a ,Wp]z — for all index pairs a, (3 — r,i,j,k. 
By bilincarity we then have [U, V]z — for any two real symmetric matrices U, V in the linear span C 
of I3, W r , Wi, Wj, W k . Note that these 5 matrices are linearly independent, because B is not reducible. 
Define the linear subspace Sw = {U € C\Uz = 0} and let S z C R 3 be the linear span of the vectors 
z = I3Z, W r z, WiZ, WjZ, WkZ. Obviously we have 

dim6>iy + dim5 2 = 5. (11) 

For any matrices U £ Sw and V G C we have U(Vz) = V(Uz) = 0. It follows that for any matrix 
U € Sw and any vector y € S z we have Uy — 0. By the same reasoning as in the previous section we 
can lead the cases dim 6^ = 3 and dimSj = 2 to a contradiction with (fTTjl . Therefore dim<S z = 1, and 
z is a common eigenvector of all matrices in C. This yields two independent linear conditions on the 
matrices in C, but C has only codimension 1 in 5(3). Therefore the considered case is not possible. 

1.2. The components of Z span the whole space M 3 . 

We have Y*Y = WW* - W*W + ZZ*. It follows that Re(Y*Y) = Re(ZZ*) and Im(Y*Y) = 
— Im(W*W). The coefficient c\ of the characteristic polynomial of Y*Y is therefore a homogeneous 
polynomial of degree 4 in the components of the elements of W, Z. By Lemma 15.21 the parts that 
depend on W and on Z are separated. It is not hard to verify that Z enters only in the form of squares 
of the elements of the imaginary part of ZZ*. If we replace these elements by corresponding terms 
depending on W using the relation Im{ZZ*) = — Im(WW*), then the polynomial simplifies to zero. 
Hence the derivative of the characteristic polynomial of Y*Y at zero is zero. Since Y*Y is PSD, and 
its characteristic polynomial has three real nonnegative roots, it must have a multiple eigenvalue at 
zero, and hence its rank cannot be bigger than 1. Therefore both Y and Z are vectors. 

Since the real parts of the products WW and WW are equal, the real parts of ZZ* and Y*Y must 
also be equal. But ReY*Y = ReY T Y, and hence by Lemma 15.11 there exist unit norm quaternions 
( z ,(y such that ZQ = C, y Y T = (( y Y) T . Since multiplication of Z from the right and Y from the 
left by unit norm quaternions does not change the matrices Z Z* ,Y*Y , we can assume without loss of 
generality that Z = Y T . 

Inserting this in (fTU|) and combining with ©, we get [W a , Wp] + Z a Zj — Zi3Z'^ = for all index pairs 
a, (3 = r, i, j, k. Let us define a linear map Wz from the linear span C of the matrices W r , Wi,Wj,Wk to 
K 3 by W z {W a ) = Z a ,a = r,i,j,k. By bilinearity we then have the relation [U, V] + Wz(U)Wz(V) T ~ 
Wz(V)Wz(U) T = for all matrices U, V <E C. Note that W z has a 1-dimensional kernel. Let U € C 
be a generator of this kernel. We have Wz(U) = and hence [U,V] = for all V € C. Therefore 
U commutes with all matrices in the 5-dimensional subspace of 5(3) which is generated by C and I3. 
But then U has to be proportional to -Z3, which contradicts the non-reducibility of B. Hence this case 
cannot occur. 

Similarly to case 1, we can treat the case when Y*Y has rank at most 1. 
2. Both Y*Y and ZZ* are of rank 2. 

Let v z , v y be generators of the kernels of ZZ* , Y*Y, i.e. vectors that satisfy v*ZZ* = 0, Y*Yv y = 
0. Consider the equation v'*(WM — MW)v y = 0, where M is a real symmetric matrix. One sees 
immediately that M — I 3 satisfies this equation. Furthermore, by ([9]), (fT0|) and Re(ZZ*) = Re(Y*Y) 
we have 

2v*{WW r - W r W)v y = 2v* z {i[Wi,W r } + j[W 3 ,W r ] + k[W k ,W r ])v y = v*{-ZZ* +Y*Y)v y = 0, 
2v* z {WW l - WiW)v y = 2v*{-[Wi,W r ] +j[Wj,Wi] - k[W t ,W k ])v y = v* z (-ZZ*i + iY*Y)v y = 0, 
2v* z {WWj - WjW)v y = 2v* z {-[Wj, W r ] - i[Wj, Wi] + k[W k , Wj])v v = v* z (-ZZ*j + jY*Y)v y = 0, 
2v* z (WW k - W k W)v y = 2v* z (-[W k , W r ] + i[W i: W k ] - j[W k , W 3 ])v y = v* z (-ZZ*k + kY*Y)v y = 0. 

Hence there exists a 5-dimensional linear subspace in 5(3) that satisfies v*(WM — MW)v y = 0. This 
subspace contains the identity matrix and thus intersects the PSD cone S+(3). But then it contains a 
rank 1 PSD matrix ££ T by Dines theorem [4]. Here £ € K 3 is a nonzero vector. 

2.1. vu=ev y =o. 

-2 



Then there exist nonzero vectors £ y , ( z such that ZC, Z = £ and Y*£ y — £■ Let a — min(|C z | , \( yl 

CyCy j 



Then we have Z = ZZ* - a^ T = Z(I 2 - a( z (*)Z* y 0, Y = Y*Y - a^ T = Y*(I 2 - a( y Q)Y y 0. It 



follows that 

R-f° °^l = ^ /3 W loR-(° \-( 1 - W 

1 a^ T J W WW + Z J ~ ' I a^ T J I W WW + Y 
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The matrix which is subtracted from B and B is clearly K 6 <g> S' + (3)-separable and hence represents 
an element of £6,3- The difference is again a PPT matrix and hence represents an element in Tq^. 

We have represented the element B as convex combination of a nonzero extremal element of £5,3 
and another element of ^,3. Since B is extremal in r 6j3 , it must be contained in £5,3- 

2.2. v* z £ ^ or £ T v y + 0. 

We have v* z (W^ T -^ T W)v y =0. Without loss of generality, let v* z £ ^ q Define q = {v* z ^)- 1 v* z W^. 
(In the case £, T v y ^ we define q — £, T Wv y (t; T Vy)^ 1 and proceed similarly.) Inserting {v* z £)q for v* z W^, 



we obtain £ J (W — ql; 



s)v y = 0. Hence we find a vector ( y € W such that £ (W — ql 3 ) + (*Y = 0, or 
(W — <?/ 3 )£ + Y*(y ~ 0- Moreover, £ commutes with quaternion scalars and we have v*(W — ql^)^ = 
v* z W£ - vlivlCl^iKWOt, = 0- Hence we find a vector ( z e H 2 such that (W - ql 3 )£ + Z( z = 0. Let 
a = min((|£| 2 + |Cz| 2 ) -1 , (|£| 2 + I Cj, 1 2 ) ^ 1 ) - Then we have 



B -a 



h 
w z 



h 
W Z 



h 0, 



B - a 



(it 



h 
w 





Y* 



(y 



Cy 



h 
w 





y * 



>- 0. 



Thus, as in the previous case, we have represented B as convex combination of an extremal element 
of ^6,3 and a nonzero element of r 6j3 . By the extremality of B in r 6j 3 we obtain B e £6,3- 

We have proven that any extremal ray of T6,3 is in £6,3- We obtain the following theorem. 

Theorem 5.5. r 6;3 = £ 6j3 . □ 
Thus a matrix in Q(2)®<S(3) is <5+(2)(g)S'+(3)-separable if and only if it satisfies the PPT condition. 
However, positivity alone is not sufficient for separability, as the following example shows. The matrix 
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G C(2)®5(3) 



is PSD, but its transpose and hence its partial transpose is not. 

Corollary 5.6. Let m > 2, n > 3. Then the cone of PSD matrices in Q(m) ® S(n) is strictly larger 
than the cone of PPT matrices. □ 



6 Counterexample against T4 4 = E 



4,4 



In this section we present an element of that is not in £4,4. 

Recall that 1^^ is the cone of blockwise symmetric complex hermitian In x 2n PSD matrices, while 
£4^ is the cone of H + (2) £g> S+ (n)-separable matrices. 

Consider the matrix B = VV* € H+(8) with 



V 



h 
w z 



w 



It is not hard to check that 
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W T , 



( 1 \ 

3i 


V / 



is symmetric and hence B is blockwise symmetric. 

Let us show that B is not H+(2) ® S'+(4)-separable. It suffices to show that there does not exist a 
separable rank 1 matrix 5 such that B — S >z 0, i.e. there is no nonzero vector v € C 5 such that Vv is 
representable as a single tensor product in C 2 <g> M 4 . 

Let us assume such a vector v exists. If its first 4 entries are zero, then Vv is not separable because 
Z is not a multiple of a real vector. Hence the first 4 entries of v form a multiple of a nonzero real 
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vector v\ and the last entry is a complex number v 2 = v 2 + iv 2 . Without restriction of generality we 
can assume that the first 4 entries of v actually form V\ . We have Wv\ + Zv 2 = Xvi for some complex 
number A = A 1 ' + i\ l . Resolving with respect to vi,v 2 ,v 2 and separating the real and the imaginary 
part, we obtain the homogeneous linear system 



ReW - \ r h ReZ -ImZ 
ImW - X l h ImZ ReZ 



Vl 



Hence the 8x6 coefficient matrix on the left-hand side must be rank deficient. The determinant of the 
matrix formed of the last six rows is proportional to A r (2 — A r ), hence A r = or A r = 2. In the first 
case adding and in the second case subtracting rows 2 and 3 leads to v 2 — 0. This yields the reduced 
systems 
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which are easily seen to have no nontrivial solution for any A 1 . 
Hence B is not separable and we obtain the following theorem. 



□ 



Theorem 6.1. T 4 4 ^ £ 4 4 . 

Theorem 16.11 and Lemma 13.11 yield the relation T 4i „ ^ S 4j „ for any n > 4 and thus the following 
result. 

Corollary 6.2. The positivity condition in the space TL{2) ® S(n) is not sufficient for H + {2) <S> S+(n)- 
separability for n > 4. □ 

Remark: Woronowicz already presented in [T5] a PPT matrix in H + (8) which is not H + (2)(g)H + (4)- 
separable. However, this matrix was not blockwise symmetric and hence not a counterexample against 
the equality T 4 . 4 = S 4i4 . 

Theorem 16.11 and Lemma 13-11 yield the relation r6, n 7^ £<6,n for any n > 4 and thus the following 
result. 

Corollary 6.3. The PPT condition in the space Q(2) ® S(n) is not sufficient for Q+(2) (g> 
separability for n > 4. □ 



7 Positivity implies *S' + (2) <g> Q+(^)-separability 

In this section we investigate the spaces 5(2) <g) Q(n) for n > 1. The dimensionality argument Lemma 
12.61 will be essentially sufficient to prove that the positivity of a matrix is equivalent to S+ (2) ® Q + (n)- 
separability. 

Let P n be the intersection of the space 5(2) Cg) Q(n) with the cone Q+(2n) and let £„ be the 
5+ (2) ® Q + (n)-separable cone. By Proposition 12.91 we have £„ C P n for any n > 1 and trivially 
Pi = Ei. Any matrix B G S(2) ® Q(n) can be partitioned as 

»-(«;*:)■ « i2 » 

where P>n, B i2 , B 2 2 € Q(n). 

Lemma 7.1. Lei B generate an extreme ray of P n . Then the rank o/B is a£ most n. 
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Proof. We apply Lemma \2~M With the notations of this lemma, define K to be the cone Q+(2n) and 
L the space S(2) ® Q(n). Since B generates an extreme ray of P n , the dimension of its face in Q + (2n) 
can be at most (8n 2 — 2n) — (6n 2 — 3n) + 1 = 2n 2 + n + 2. Now suppose that B has rank at least 
n+ 1. Then the dimension of its face in Q + {2n) equals at least the dimension of Q(n + 1), which is 
2{n + l) 2 — (ra+ 1) = 2ri 2 + 3n + 1 > 2n 2 + n + 2. This contradiction completes the proof. □ 

Let us now define a group homomorphism H® : GL n (W) -> Aut(Q+(n)). It shall assign to any 
regular quaternionic n x n matrix S the automorphism A t— > SAS*. Note that for any S € GL„(H), 
the map id 5 ( 2 ) ®H^(S) is both in Aut(P n ) and in Aut(S n ). 

Lemma 7.2. Let A, B 6 Q + (n). TTien i/iere exists a matrix S G GL n (H) smc/i £/iai both SAS* and 
SBS* are diagonal. 

Proof. Assume the conditions of the lemma. Let I < n be the rank of A + B. Since A + B >z 0, there 
exists a regular matrix S% such that Si(A + B)S* is equal to diag(/;, 0„_;). Now note that SiAS* ^ 
Si(A + B)Sl, SxBSl < Si(A + B)Sf and hence SxAS^ = diag(Ai, n _,), SxBS* = diag(Si,O n _0 
for some matrices A\,Bx £ Q + (l). Let now U be a hyperunitary I x I matrix such that U(A\ — 
Bi)U* is diagonal (see Proposition IA.2I in the appendix). Then one easily sees that the matrix S = 
diag(f/, J„_;)5i satisfies the assertion of the lemma. □ 

Corollary 7.3. Let BeP„. Then there exists S € GL„(IHI) such that both the upper left and the lower 
right n x n block of the matrix B' = (ids(2) (8>i?^(5))(B) are diagonal. □ 

Lemma 7.4. Let B generate an extreme ray of P n and be partitioned as in ilty) . Then there exists 
S E GL n (M) such that the map H^(S) diagonalizes all three blocks 5n, -B12, i?22- 

Proof. We prove the lemma by induction. For n = \ the assertion of the lemma holds trivially. Let us 
now assume that it holds for n — 1. 

Assume the conditions of the lemma. By the preceding corollary we can assume without restriction 
of generality that B\\ and B22 are already diagonal. 

Suppose that the rank of B21 is strictly smaller than n. Then at least one diagonal element of B22, 
say the last one, is zero. Since B ^ 0, the last row and the last column of B are zero. But then also the 
last row and the last column of Byi are zero. Hence all three matrices Bn, B12, B22 are block-diagonal, 
with an upper left block of size (n — 1) x (n — 1) and a lower right block of size lxl. By assumption of 
the induction there exists 5" £ GL n _i(M) such that the map /f^L 1 (5) diagonalizes all the upper left 
blocks. It then follows that the matrix S = diag(S", 1) satisfies the assertion of the lemma. 

Let us now assume that B22 has full rank. Then there exists Si € GL n (M) such that S1B22S* = L n . 
Define B' = (id 5(2) (8)F"(5i))(B) and B' n = SxB n S$, B[ 2 = SiB 12 Si*. Let now U be a hyperunitary 
n x n matrix such that UB' 12 U* is diagonal. By Lemma [7.11 the rank of B' equals n. This implies 
B' u = (B[ 2 ) 2 - It follows that UB'^U* = {UB' 12 U*) 2 is also diagonal. Thus the matrix S = USi 



satisfies the assertion of the lemma. □ 

Corollary 7.5. Let B generate an extreme ray of P n . Then B 6 £„. □ 
This yields the following result. 

Theorem 7.6. P n = £„ for any n > 1. □ 



8 Conclusions 

Let us summarize our results. As mentioned in the introduction, the cone of S+(m) (E> 5+(n)-separable 
matrices in 5(m)(8<S(n) coincides with the cone of PSD matrices in this space if and only if min(n, m) < 
2. This is a consequence of the results in [TD]. The cone of H + (m) ® H + (n)-separable matrices in 
H(mn) = W(m)(g>H(n) coincides with the cone of PPT matrices in Tt(mn) if and only if min(n, m) = 1 
or m + n < 5. This is a consequence of the results in [T5]. We illustrate this in the following tables. 



m\n 


2 


3 


> 4 


m\n 


2 


3 


> 4 


2 


PSD 


PSD 


PSD 


2 


PPT 


PPT 


N 


3 


PSD 


N 


N 


3 


PPT 


N 


N 


> 4 


PSD 


N 


N 


> 4 


N 


N 


N 



Semidefinite descriptions of the Semidefinite descriptions of the 

iS+(m) (g) S+ (n)-separable cone H+{m) ® i?+(n)-separable cone 

Here "PSD" indicates that the cone of separable elements in the space corresponding to the pair 
(to, n) equals the cone of positive semidefinite matrices, and "PPT" indicates that the cone of separable 
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elements equals the cone of PPT matrices. An "N" indicates that the separable cone is described by 
neither the PSD cone nor the PPT cone. 

Theorems 13.91 l3~12l Corollary 16.21 and Lemma l3~2l lead to the following result. 

Theorem 8.1. The cone of H + (m) (g) S+(n)- separable matrices in H.(m) (g> 5(n) coincides with the 
cone of PSD matrices in this space if and only if n < 2 or m — 1 or m + n < 5 . □ 



Theorem 17. 61 Lemma T3. 21 and Corollary [576] lead to the following result. 

Theorem 8.2. The cone of Q+(m) ® S+(n)- separable matrices in Q(m) (g> S(n) coincides with the 
cone of PSD matrices in this space if and only if n < 2 or m = 1. □ 



Theorem l5.51 Corollary 16.31 and Lemma 13721 lead to the following result. 

Theorem 8.3. The cone of Q+(m) (g> S+(n)- separable matrices in Q(m) <g> S(n) coincides with the 
cone of PPT matrices in this space if and only if m < 2 and m + n < 5. □ 
We can summarize these results in the following tables. 



m\n 


2 


3 


> 4 


m\n 


2 


3 


> 4 


2 


PSD 


PSD 


N 


2 


PSD 


PPT 


N 


3 


PSD 


N 


N 


3 


PSD 


N 


N 


> 4 


PSD 


N 


N 


> 4 


PSD 


N 


N 



Scmidcfinitc descriptions of the 
H + (m) ® S+ (n)-separable cone 



Scmidefmite descriptions of the 
Q + (m) (g> 5+(n)-separable cone 



In all four tables, for min(n, m) = 1 the cone of separable matrices is trivially equal to the cone of 
PSD matrices. 

The preceding two theorems yield also the following nontrivial result. 

Corollary 8.4. A matrix in Q(2) ® 5(2) is PSD if and only if its transpose is PSD. □ 
Apart from the trivial cases Q(2) = Q(2) <g> 5(1) and S(n) = Q(l) ® S(ri) the space Q(2) ® 5(2) is 
thus the only tensor product space Q(m) <g> 5(n), m,n6 N+, where the positivity property is invariant 
with respect to transposition (cf. Corollaries IA.8I and IA.9I in the appendix) . 



References 

[1] Hel Braun and Max Koecher. Jordan-Algebren, volume 128 of A Series of Comprehensive Studies 
in Mathematics. Springer, Berlin, New York, 1966. 



[8 
[9 

no 



J.L. Brenner. Matrices of quaternions. Pacific J. Math., 1:329-335, 1951. 

M.-D. Choi. A Schwarz inequality for positive linear maps on C-algebras. Illinois J. Math., 
18(4):565-574, 1974. 

Lloyd L. Dines. On linear combinations of quadratic forms. Bull. Amer. Math. Soc, 49:388-393, 
1943. 

Leonid Gurvits and Howard Barnum. Separable balls around the maximally mixed multipartite 
quantum states. Phys. Rev., 68:042312, 2003. 

Ryszard Horodecki, Pawel Horodecki, and Michal Horodecki. Separability of mixed states: neces- 
sary and sufficient conditions. Physics Letters, 223(l):l-8, 1996. 

H.C. Lee. Eigenvalues and canonical forms of matrices with quaternion coefficients. Proceedings 
of the Royal Irish Academy, Sect. A, 52:253-260, 1949. 

A. Peres. Separability criterion for density matrices. Phys. Rev. Lett., 77:1413-1415, 1996. 

Erling St0rmer. Positive linear maps of operator algebras. Acta Mathematica, 110:233-278, 1963. 

F.J. Terpstra. Die Darstellung biquadratischer Formen als Summen von Quadraten mit Anwen- 
dung auf die Variationsrechnung. Mathematische Annalen, 116:166-180, 1939. 



[11] E. B. Vinberg. The theory of convex homogeneous cones. Transactions of Moscow Mathematical 
Society, 12:340-403, 1963. 



19 



[12] Frank Warner. Foundations of Differentiable Manifolds and Lie Groups. Springer, New York, 
Berlin, Heidelberg, Tokyo, 1983. 

[13] N.A. Wiegmann. Some theorems on matrices with real quaternion elements. Can. J. Math., 
7:191-201, 1955. 

[14] L.A. Wolf. Similarity of matrices in which the elements are real quaternions. Bulletin of the 
American Mathematical Society, 42:737-743, 1936. 

[15] S.L. Woronowicz. Positive maps of low dimensional matrix algebras. Reports on Mathematical 
Physics, 10:165-183, 1976. 

[16] V.A. Yakubovitch. Factorization of symmetric matrix polynomials. Doklady Akademii Nauk SSSR, 
194(3):1261-1264, 1970. 

[17] Fuzhen Zhang. Quaternions and matrices of quaternions. Linear algebra and its applications, 
251:21-57, 1997. 



A Quaternions 

In this section we provide some basic facts about quaternions and quaternionic matrices, which will be 
used in the paper. Most concepts from the theory of real and complex hermitian matrices carry over to 
the quaternionic case. For a detailed treatment of the properties of quaternionic matrices see [2"] .[17). 

We denote the algebra of quaternions by H, with the usual product rules 

i 2 = j 2 = k 2 =ijk= -1 (13) 

for the generators. This algebra is not commutative. However, any non-zero quaternion q = q r + iqi + 
jqj + kqk has a two-sided inverse q^ 1 = q/\q\ 2 , where q = q r — iqi — jqj ~~ kqk is the quaternionic 

conjugate of q and \q\ = ^Jq 2 + qf + q 2 + q 2 . is the norm of q. The real numbers q r ,qi,qj,qk are the 

components of q. The norm has the properties \q\ = \q\ and \pq\ = \p\\q\- The conjugate of a product 
is given by pq = qp. Denote by Re q the real part (q + q)/2 of q. 

Let Qij : H — > H be the map defined by Qij : q = q r + iqi + jqj + kqk i— > q r ~ iqi ~ jqj + kqk- 
It is not hard to see that Qij is an automorphism of the quaternion algebra, since it leaves relations 
(I13p invariant. Likewise we can define the automorphisms Qik and Qjk, which switch the signs of the 
corresponding components. There exist also other automorphisms, e.g. the cyclic permutation of the 
imaginary units i, j, k. 

Let us now consider vectors and matrices of quaternions. 

Denote the space of column vectors with n quaternionic entries by H™ and the space of quaternionic 
matrices of size m x n by W mxn . The norm ||u|| of a quaternionic vector v can be defined as y/v*v. 
Here v* is the quaternionic conjugate transpose of v. This norm equals the Euclidean norm in the 
associated real vector space R 4 ™ of the components of the elements of v. 

The conjugate A of a quaternionic matrix A is defined entrywise. Let A* denote the conjugate 
transpose of A. Then for matrix products we obviously have (AB)* — B* A* . 

The quaternionic matrix algebras have a representation as real matrix algebras. Namely, if we write 
quaternionic matrices A as A = A r + iAi + jAj + kAk, where A r , Ai, Aj, A^ are real matrices, then the 
isomorphism is given by 

/ A r —Ai ~Aj —Ak \ 
Ai A r Ak Aj 
Aj Aj~ A r Ai 
\ A/~ Aj Ai A r J 

The conjugate transpose of a quaternionic matrix corresponds to the transpose of its real counterpart, 
l q (A*) = {l q {A)) T - 

The trace of matrix products is not invariant under cyclic permutations, even for scalars. Never- 
theless, we have the following weaker assertion. 

Proposition A.l. For quaternionic matrices A, B of appropriate sizes the relation Re(tr(AB)) = 
Re(tr(BA)) holds. 

Proof. For a square matrix A we have trl q (A) = ARe{trA). Hence Re(tr(AB)) = \tr{l q (A)l q (B)) = 
\tr{l q {B)l q {A)) = Re{tr{BA)). □ 



1 q : A 



(14) 
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Let A be an n x n square quaternionic matrix. If there does not exist a nonzero quaternionic vector 
v such that Av = 0, then there exists an inverse quaternionic matrix A -1 with AA^ 1 = A~ 1 A = I n 
[5] • We will call such matrices regular. A regular matrix A can be viewed as an R-linear automorphism 
of H", which maps a vector v to the vector Av. If there exists n ^ such that Av = then we 
call A singular. Obviously I q (A^ 1 ) = (T q (A))~ 1 , and the inverse A^ 1 of a regular matrix is unique. 
The inverse of a product is given by (AB)^ 1 = B~ 1 A~ 1 , the inverse of the conjugate transpose by 
(A*) -1 = (A^ 1 )* , which is checked easily by applying the isomorphism T q . 

If for a square quaternionic matrix A and a quaternionic vector v there exists a quaternion A such 
that Av — vX, then we call v right eigenvector of A and A the corresponding right eigenvalue. An nx n 
quaternionic matrix has n right eigenvalues, which are determined up to similarity transformations 
(i.e. transformations A i— > qXq^ 1 for q E H, q ^ 0). The right eigenvalues are preserved under similarity 
transformations of the matrix (i.e. A i— » SAS^ 1 for regular S). HA is triangular, then the diagonal 
elements are representatives of the (left and right) eigenvalues [7]. 

If for a square matrix U we have UU* = I, then we call U hyperunitary. A matrix U is hyperunitary 
if and only if the matrix T q (U) is orthogonal, hence the hyperunitary matrices of size n x n form a 
compact group, the compact symplectic group Sp(n). Hyperunitary matrices can be viewed as norm- 
preserving M-linear automorphisms of i™ (though not every such automorphism can be represented 
by a hyperunitary matrix). 

Any quaternionic matrix A can be decomposed as A = UDV, where U, V are hyperunitary matrices 
and D is a diagonal matrix, which has the same size as A, with real nonnegative entries. The diagonal 
elements are called singular values of A [7], [13] . For a quaternionic matrix A, we call the number of 
its positive singular values the rank of A. Obviously an n x n matrix is regular if and only if it has 
rank n. The rank is invariant under multiplications by regular matrices from the left and from the 
right [IT] . It can be shown that the rank of a quaternionic matrix is equal to the maximal number of 
right linearly (over H) independent columns and left linearly independent rows of the matrix [14j . 

We call a square quaternionic matrix A hermitian if A = A*. The hermitian matrices of size 
n x n form a (2n 2 — rt)-dimensional vector space Q(n) over the reals. The real part A r of a hermitian 
quaternionic matrix is symmetric, while the three imaginary parts Ai,Aj,Ak are skew-symmetric. 
There is a scalar product on Q(n) given by 

(A, B) = tr(A r B r - AiBi - AjBj - A k B k ) = Re(tr(AB)) = Re(tr(BA)). 

Proposition A. 2. Let A E Q(n). Then there exists U E Sp(n) such that A — UDU* , where D is 
a diagonal matrix with real entries. The columns of U are right eigenvectors of A, while the diagonal 
elements of D are the corresponding right eigenvalues. 

This follows from a well-known result on factorization of quaternionic matrices [2] , namely that for 
any square quaternionic matrix A there exists a hyperunitary matrix U such that A — UTU*, where 
T is upper triangular. Namely, if A is hermitian, then T must also be hermitian and hence diagonal 
and real. 

The right eigenvalues of a hermitian matrix are real and uniquely determined. The singular values of 
a hermitian matrix are given by the absolute values of its eigenvalues, hence its rank equals the number 
of non-zero eigenvalues. Since the eigenvalues A; of A are well-defined, its characteristic polynomial 
p(X) = nlLi(^ ~ ^i) is well-defined and has real coefficients. 

Proposition A. 3. Let A E Q(n). Then the characteristic polynomial ofT q (A) is the fourth power of 
the characteristic polynomial of A. 

Proof. Consider the matrix T q {A) as an element of Q(4n). Its eigenvalues are real and do not change 
under conjugation with the hyperunitary matrix 





/ 1 i j k 


\ 


1 


—1 i —jk 


2 


j k 






\ 1 -i -j k 


) 



This conjugation leads to a block-diagonal matrix with blocks A, A r + iAi — jAj — kA k ,A r — iAi + 
jAj — kAk,A r — iAi — jAj + kA k . Hence the characteristic polynomial of T q {A) equals the product 
of the characteristic polynomials of these 4 blocks. But the last three blocks are the images of A 
under the automorphisms Qij, Qjk, Qik- Since the real line is invariant under these automorphisms, the 
eigenvalues of these blocks coincide with the eigenvalues of A. Therefore their characteristic polynomials 
also coincide with that of A. This completes the proof. □ 

It follows that the characteristic polynomial of a hermitian matrix A can be computed by taking 
the unique fourth root of the characteristic polynomial of T q (A) whose leading coefficient equals 1. 
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A matrix A S Q(n) is called positive semidefinite (PSD) if for any vector v £ HP we have v* Av > 
and positive definite (PD) if for any non-zero vector v £ HP 1 we have v* Av > 0. Note that (v*Av)* = 
v*A*v, hence v* Av £ R. Clearly A is PSD if and only if all of its eigenvalues are nonnegative and PD 
if and only if all of its eigenvalues are positive. Let us denote the cone of PSD matrices in Q{n) by 
Q+(n). 

For any A £ Q+(n) of rank I there exists V £ EP xi of rank I such that A = VV* . Namely, let 
A = U ( - jj J U*, where U is hyperunitary and D is a real diagonal I x / matrix with positive 



£1/2 



diagonal elements. Then we can choose V — U y ^ j . We have the freedom of multiplying the 

factor V by a hyperunitary matrix from the right. On the other hand, if A = VV* , then A y and the 
rank of A equals the rank of V. Let us define A 1 / 2 = UD 1 / 2 U*, then A 1 / 2 £ Q+(n) and A 1 ! 2 A 1 ! 2 = A. 
If A is PD, then A 1 ! 2 is also PD and its inverse is given by yJ~ 1/2 = [7D~ 1/2 [/*. 

If for some vector v £ H™ and a matrix A £ Q+(n) we have v* Av — 0, then Av = 0. This is because 
v*Av = (A 1 / 2 v)*(A 1 / 2 v) = l^ 1 / 2 !-! 2 = yields A 1 / 2 v = and therefore A 1 / 2 A 1 / 2 v = 0. 

Now let n — n\ + ri2- If A = ( ^f 11 ] £ Q + (n) is partitioned in four blocks such that 

V A21 A22 j 

An £ Q{ni) and A22 £ Q{n2) are square matrices, and v £ Ti ni is such that Anv — 0, then we have 

0n 2 X 1 



also A21V — 0. This can be seen as follows. If we define v — [ _ I , then v* Av = v*Auv = 



Hence Av = ( J \ llV I = 



A21V 

Let a C {1, . . . , n} be a subset of indices, let A = VV* be a PSD nxn matrix. Then the submatrix 
A of A that consists of those elements whose row and column numbers are in a is also PSD, because it 
can be represented as product VV* , where V consists of those elements of V whose row numbers are 
in a. Moreover, if A y 0, then A y 0. On the other hand, a block-diagonal hermitian matrix is PSD 
if and only if each block is PSD, and it is PD if and only if each block is PD. 

It is well-known that the cone Q+(n) is convex and self-dual [llj . Any automorphism of the 
quaternion algebra induces an automorphism of the cone Q+(n), because it leaves the real line invariant. 

The cone i?+(n) of complex hermitian PSD matrices is an intersection of Q+(n) with a linear 
subspace. Namely, if A is hermitian with zero j- and fc-components, then A £ Q+{n) if and only if 
1q(A) £ S+(4n) and hence if and only if A £ H + (n). 

Any PSD quaternionic matrix can be written as sum of matrices of the form vv* , where v are 
vectors. These vectors can be obtained e.g. as columns of the factor V of A = VV* . Thus Q+(n) is the 
convex conic hull of the set of rank 1 matrices {vv* \ v £ H™}, and any such rank 1 matrix generates 
an extreme ray of Q + (n). The matrix vv* , where v £ H™ is a nonzero vector, has only one positive 
eigenvalue, namely \v\ 2 = v*v. Hence I — vv* y for all v in the unit ball. 

Lemma A. 4. Let S £ H nx ' be of rank I, where I < n, and let A £ M lxL . Then SAS* = is equivalent 
to A = 0. If A £ Q(l), then SAS* h is equivalent toAyO. 

Proof. Clearly A = implies SAS* = and Ay{) implies SAS* b 0. Let S = UDV be the singular 
value decomposition of S. Then S* S = V* D T DV . But D T D is PD, hence S*S is also PD and 
invertible. Then we have A = (S*Sy 1 S*(SAS*)S(S*S)~ 1 . Hence the converse implications hold 
too. □ 

Lemma A. 5. Let S £ M nxl be of rank r. Then there exists T £ M^ n ~ r " >xn of rank n — r such that 

TS = 0. Moreover, if for some S' £ W nxl ' we have TS 1 = 0, then there exists W £ W lxV such that 
S' = SW. 

Proof. Let S = U ( „ „ ® rx( . l ~ r ) ) 1/ be the singular value decomposition of S, where D 

\ "(n— r)xr U(„_ r ) x (;_ r ) J 

is a positive definite diagonal r x r matrix. Define T = (0(„_ r ) Xr I n - r )U* . Then T is of rank n — r 
and TS = 0. Suppose S' satisfies the assumptions of the lemma. Then the last n — r rows of U*S' are 
zero. Denote the matrix given by the first r rows by P. Define 



W = V* 



D- X P 

0(l-r)xl' 



Then we have 



D Orx(l-r) \ ( D-^P \ TT ( P 



SW = U[ - "rxy-r; - =f/ U )=U{U*S') = S'. □ 
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Let us investigate the facial structure of Q+(n). 

Proposition A.6. Let S G W ixl be of rank 1,1 < n. Then the set {SBS* \ B > 0} is a face ofQ + (n). 
This face is isomorphic to Q + (l). Any face of Q + (n) can be expressed in such a way. 

Proof. Let A G Q+(n) be arbitrary. Denote the rank of A by I. Then there exists S G M nxl of rank 
I such that A = SS*. By Lemma [PI the map B i— ► SBS* is an isomorphism between the space 
and the linear subspace L = {SBS* \B G 2(0} C Q(^), and maps to the intersection 

L + = LP Q+(n). Moreover, A is in the interior of L + . Let us show that L + is a face of Q+(n). 

Let Mi,M 2 G Q+(n) such that Mi±Mz = SBS* for some B G Q+(Z). We have to show that 
M U M 2 G L+. Define C = M 2 - M x G Q(n), then S B S* +aC £ Q+ (n) for all a G [-1/2, +1/2]. By the 
previous lemma there exists T G H^)*™ of rank n-l such that T5 = 0. We have T(SBS*+aC)T* = 
aTCT* h for all a G [-1/2, +1/2]. Hence TCT* = and T{SBS* + aC)T* = for all a G K. 
But ,555* + «C h for a G [-1/2, +1/2], therefore we also have (SBS* + aC)T* = aCT* = and 
TC — 0. By the previous lemma there exists a matrix W such that C = SW. Since C is hermitian, 
we also have C = W*S*. It follows that W = {S*S)~ 1 S*SW = (S* S)~ 1 S*W* S* and therefore 

C* = S(S*S)~ 1 S*W*S* = SWS (S* S) - 1 S* = s {s ' srls * w "+ ws< > s * srl S*. We have constructed a 
hermitian matrix B' such that C = SB'S* G L. Therefore M 1 ,M 2 € L and hence M 1 ,M 2 G L+. 

Let now S G HP xi be given and of rank /. Then we define A = SS* and proceed as above. □ 

Similar results on the facial structure of S+(n) and H+(n) can be obtained by the same line of 
reasoning. 

Lemma A. 7. Let n = n% + n 2 . Let A = ( ^ n ^f 12 ) G QW fre partitioned in four blocks such 

V ^21 ^22 / 

</ia< An G is PP. T/ien A ^ if and onZy if A 22 — A 2 \A~^ A\ 2 y 0. Moreover, the rank of A 

is the sum o/rk(A22 — A 2 \A^ A\ 2 ) and rk An. 

Proof. If A 22 - A 21 A^A 12 is PSD of rank I, then there exists Z G H" 2X| of rank I such that 

A= ( Ai 2 o ni xi\( < 2 o nixl y 

^ A 21 A-, 1/2 Z ) \ A 21 A U 1/2 Z J 

Hence A is PSD. 

Let A y and let x 2 G H" 2 be arbitrary. Define £1 = —A~ 11 A\ 2 x 2 G H™ 1 , then we have 

(a;* a^)^ f ^ J = x 2^2iA^A 12 x 2 - x^A 21 A^A 12 x 2 - x^A 21 A^A 12 x 2 + £2^22^2 
= £2(^22 - AaiA^Aia)^ > 0, 

Now let A}f = J/PP* and Z = U'W be the singular value decompositions of A\{ 2 and Z. 



and A 22 - ^i^n 1 ^ t 

No- 
Define 



tf* \ { A\[ 2 niX i\(U \_( D nixl x x(ni+0 



and partition it into matrices Pi G H( ,ll +') x (" 1+i ), P 2 G H^-^O^+O. Note that P and Pi are lower 
triangular matrices, with positive elements on the diagonal. Hence P* is hence invertible. But then 
P*Pi and P* P = P*P% + P 2 P 2 are PD. Hence the rank of P equals the number m + 1 of its columns, 
which in turn is equal to rk An + rk(A_22 — A 2 iA^Ai 2 ). This completes the proof. □ 

Corollary A. 8. A matrix A = (a a p) G Q(2) is P5T> if and only if an > 0,a 22 > 0, a\\a 22 > |a.i2 1 2 - 
In particular, the cone Q+(2) is invariant with respect to transposition. □ 
However, in general A y does not yield A y 0. For example, we have 



yQ, i i 1 fc i = -3 < 0. 





Corollary A. 9. P/ie cone Q+(n) is invariant with respect to transposition if and only if n < 2. □ 
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